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Abstract 


A  displacement  controlled  dynamic  test  system  and 
instrumentation  for  dynamic  fracture  testing  of  concrete 
specimens,  was  developed.  Also  a  hybrid  experimental-numerical 
procedure  for  evaluating  the  dynamic  fracture  data  thus 
generated  was  developed.  The  experimental  data  consisted  of  the 
applied  load  and  transient  strains  in  single  edge  notched,  three 
point  bend  specimens  which  were  subjected  to  controlled 
displacement  loading.  In  addition,  moire  interferometry  data  was 
obtained  in  three  point  bend  specimens  subjected  drop  weight 
loading.  Also  the  applied  load,  transient  strains  and  crack 
opening  displacement  were  measured  in  crack-line  wedge-loaded, 
double  cantilever  beam  specimens  which  were  subjected  to 
controlled  wedge  displacement  loading.  These  test  systems  and 
the  hybrid  analysis  were  used  to  determine  the  strain  rate 
sensitivity  of  the  elastic  properties  and  the  fracture  process 
zone  of  rapidly  fracturing  concrete  specimens.  The  increase  in 
both  tensile  and  compressive  moduli  of  elasticity  and  the 
tensile  strength  with  increasing  strain  rate  were  quantified. 
The  fracture  process  zone  decreased  in  magnitude  with  increasing 
s^rain  rate  but  the  tensile  strength  necessary  to  propagated  the 
tip  of  the  fracture  process  zone  under  a  tensile  strength 
criterion  increased  substantially  with  the  strain  rate.  The  net 
effect  was  that  dynamic  fracture  responses  of  a  rapidly 
propagating  crack  in  concrete  resembled  that  of  linear  elasto- 
dynamic  fracture  mechanics  with  influence  of  the  fracture 
process  zone  decreasing  with  increasing  strain  rate  of  loading. 
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CHAPTER  ONE 


INTRODUCTION 


1 . 1  General  Remarks 

In  Part  I  of  this  report,  a  general  survey  of  the 
literature  on  fracture  mechanics  and  dynamic  fracture  testing  of 
concrete  was  provided.  A  constitutive  relation,  which  related 
the  crack  closure  stress  and  crack-opening  displacement  (COD)  in 
quasi-statically  fracturing  three-point-bend  and  crack-line, 
wedge-loaded,  double  cantilever  beam  (CLWL-DCB)  concrete 
specimens  was  obtained  by  a  hybrid  experimental-numerical 
procedure.  The  procedure  consisted  of  matching  the  computed  COD 
with  the  COD  measured  by  Moire  interferometry.  The  former  was 
obtained  through  an  elastic  finite  element  analysis  of  the 
concrete  fracture  specimens  with  a  postulated  constitutive 
equation  relating  the  the  crack  closure  stress  and  the  COD.  This 
constitutive  relation  was  then  used  successfully  to  model  the 
dynamic  fracture  responses  of  published  experimental  results  of 
impacted  three-point  bend  concrete  specimens. 

In  this  Part  II,  the  results  of  experiments  on  hybrid 
experimental-numerical  analyses  of  dynamically  loaded  three- 
point  bend  and  CLWL-DCB  concrete  specimens  are  reported.  For 
four  fracturing  three-point  bend  specimens,  moire  interferometry 
was  also  used  to  measure  the  transient  COD  distributions. 
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1 . 2  Objectives 


The  objective  of  this  research  is  to  delineate  the  dynamic 
fracture  mechanism  of  concrete  under  mode  I  loading.  The 
ultimate  objectives  of  this  research  are  to  obtain  an  accurate 
description  of  the  strain  rate  sensitivity  of  the  mechanical 
properties  of  concrete  and  a  numerical  simulation  of  the  dynamic 
fracture  behavior  of  concrete. 

To  achieve  those  objectives,  the  following  four  tasks  were 
completed : 

1.  A  displacement  controlled  dynamic  fracture  test  system  was 
developed . 

2.  The  strain  rate  sensitivity  of  concrete  was  established 
through  testing  two  different  types  of  specimens,  a  three- 
point  bend  specimen  and  a  crack-line  wedge-loaded,  double¬ 
cantilever  beam  (CLWL-DCB)  specimen. 

3.  A  numerical  model  of  the  propagating  crack  with  a  fracture 
process  zone  (FPZ)  was  developed. 

4 .  An  inverse  technique  was  developed  for  determining  the 
mechanical  properties  required  to  simulate  accurately  the 
global  experimental  responses. 
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1 . 3  Scop* 


Details  of  the  dynamic  loading  research  accomplished  to 
date  are  documented  in  this  report.  Those  results  are  preceded 
in  Chapter  2  by  a  literature  survey  on  dynamic  fracture  of 
concrete.  Chapter  3  describes  general  features  of  the  dynamic 
loading  system,  instrumentation  and  data  recording  systems  used 
along  with  general  details  of  the  numerical  analysis  procedures. 
Chapter  4  describes  the  findings  from  the  displacement 
controlled  three  point  bend  tests  including  details  of  the  test 
results  and  their  evaluation.  Chapter  5  describes  the  findings 
from  the  impact  loading  tests  on  three  point  bend  tests  together 
with  details  of  the  test  results  and  their  evaluation.  Chapter  6 
describes  the  results  of  the  displacement  controlled  tests  on 
CLWL-DCB  specimens.  Chapter  7  summarizes  the  investigation  and 
the  principal  conclusions  that  can  be  drawn  from  it. 
Publications  from  this  study  and  related  work  are  summarized  in 
Chapter  8  and  recommendations  for  continuation  of  this  work  are 
presented  in  Chapter  9. 
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CHAPTER  TWO 


LITERATURE  REVIEW  ON 
DYNAMICS  FRACTURE  MECHANICS  OF  CONCRETE 

2 . 1  Introduction 

Although  concrete  is  a  very  common  structural  material,  the 
state  of  science  relating  to  its  fracture,  particularly  its 
dynamic  fracture,  lags  behind  the  current  advanced  state  of 
science  for  dynamic  fracture  of  metals  and  polymers.  For 
concrete  members,  where  failure  is  controlled  by  yielding  of  the 
reinforcement  or  crushing  of  the  concrete,  common  practice 
provides  reasonably  reliable  information.  However,  for  members 
with  failure  modes  controlled  by  crack  propagation 
considerations,  i.e.,  concrete  fracture,  diagonal  tension 
failure  or  splitting  failures,  common  practice  does  not  yield 
reliable  information  since  little  is  known  about  dynamic  crack 
propagation  in  concrete  under  dynamic  loading. 

While  research  on  the  dynamic  fracture  properties  of 
concrete  is  in  its  infancy  [1,2],  the  extensive  non-linear 
effects  present  in  the  static  fracture  of  concrete  suggest  that 
the  direct  technology  transfer  of  linear  elastodynamic  fracture 
mechanics  concepts  from  brittle  materials  to  concrete  is  not 
likely  to  be  successful.  Indeed,  available  evidence  [3]  suggests 
that  dynamic  fracture  of  concrete,  in  the  absence  of  crack 
curving,  closely  mimics  dynamic  ductile  tearing  which  reduces 
the  crack  driving  force  and  hence  the  crack  velocity  in  metals 
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and  polymers  [4] .  The  similarities  between  the  apparent  dynamic 
fracture  responses  of  ductile  metals/polymers  and  concrete 
suggest  that  plastic  deformation  in  the  former  can  be  replaced 
by  aggregate  crack  bridging  in  the  latter.  The  dynamic  extension 
of  the  similarity  between  the  static  Dugdale  strip  yield  zone 
[5]  for  ductile  metals/polymers  and  the  static  fracture  process 
zone  due  to  aggregate  crack  bridging  forces  in  concrete  [6-8], 
is,  however,  not  straight  forward  due  to  the  large  increases  in 
tensile  strength  [9-12],  and  the  even  greater  relative  increases 
in  energy  absorption  [3],  reported  for  concrete  loaded 
dynamically . 

For  both  static  and  dynamic  loadings,  crack  propagation  in 
concrete  is  a  continuous  process  and  the  resulting  failure, 
unlike  fracture  in  metals  and  polymers  occurs  relatively  slowly. 
Even  for  concrete  subjected  to  static  loading,  a  crack  takes 
time  to  fully  develop  [13,14],  and  therefore  it  is  only  to  be 
expected  that,  for  dynamic  loadings,  cracks  propagate  at 
relatively  low  velocities  [15]. 

The  foregoing  observations  suggest  two  mechanistic  changes 
in  the  resistance  of  concrete  to  crack  propagation  with 
increasing  strain  rates.  One  change  reflects  the  behavioral 
change  from  a  domination  by  aggregate  crack  bridging  forces  for 
the  fully  developed  fracture  process  zone  present  in  static 
fracture,  to  a  mechanism  involving  only  a  partially  developed 
zone  for  dynamic  loading.  The  slow  crack  velocity  for  dynamic 
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loading  makes  it  distinctly  possible  that  the  energy  input  drops 
off  before  the  fracture  process  zone  is  fully  developed.  The 
second  mechanism  involves  the  development  of  relatively  more 
micro-cracking,  along  and  in  the  vicinity  of  the  line  of  the 
extending  crack,  for  dynamic  as  opposed  to  static  loading.  More 
micro-cracking  would  result  effectively  in  a  more  circuitous 
crack,  and  therefore  greater  energy  absorption  with  increasing 
strain  rates. 

2 . 2  Fracture  Behavior  of  Concrete 

Concrete  can  be  modeled  as  a  composite  material  consisting 
of  aggregates  and  hardened  cement  paste.  The  elastic  moduli  of 
the  aggregates  and  paste  are  essentially  constant  to  the  point 
of  failure  [16].  As  shown  in  Figure  2.1  the  behavior  of  concrete 
is  highly  nonlinear  due  to  bond  cracks  between  the  aggregates 
and  paste,  and  progressive  cracking  as  the  load  increases  [17] . 
The  numerous  bond  cracks  that  exist  at  the  interfaces  between 
mortar  and  coarse  aggregates  in  concrete  specimens  that  have 
never  been  loaded  are  generated  by:  water  lenses,  which  develop 
under  coarse  aggregates  during  setting  of  the  fresh  concrete; 
volume  change  in  the  cement  during  hydration;  and  dry  shrinkage 
of  hardened  concrete  [16-19].  The  stress-strain  curve  for 
concrete  loaded  in  compression  can  be  divided  into  the  four 
parts  shwon  in  Fig.  2.2  [20]  .  Within  each  part,  there  are 
different  growth  rates  for  the  bond  cracks  between  the 
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Figure 


Figure  2 . 


I 


.  1:  Stress-Strain  Diagram  of  Hardened  Cement  Paste, 
Aggregate  and  Concrete  [16]  . 


Strain 


2:  Stress-Strain  Curve  for  Concrete  under  Compressive 
Loading  [20] . 
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aggregate  and  the  matrix  and  the  cracks  within  the  matrix 
itself.  Figure  2.3  shows  cracking  maps  drawn  from  X-ray 
images  [22]  of  a  model  of  concrete  during  loading.  Those  maps 
illustrate  the  relative  developments  with  loading  of  bond  and 
matrix  cracks . 

Fracture  behavior  in  tension  is  very  important  since 
fracture  under  compression  is  actually  initiated  by  tension 
forces.  However,  only  a  few  experiments  on  the  complete  stress- 
deformation  behavior  in  tension  are  presented  in  literature  [23- 
29] .  The  early  experiments  of  Evans  and  Marathe  [24]  showed  that 
the  cracking  occurred  between  90  and  140  micro-strain,  and  that 
the  corresponding  normal  stress  is  between  68  and  89  percent  of 
the  maximum  average  stress.  These  values  were  larger  than  those 
(30  to  80  micro-strain)  reported  by  Kaplan's  [30]  because  of  the 
initiation  of  micro-cracks  within  the  gage  length.  The  strain 
reading  for  cracking  was  found  to  be  gage-length  dependent. 

Figure  2.4  shows  a  typical  example  of  the  results  from  the 
direct  tension  test  conducted  by  Reinhardt  and  his  co¬ 
workers  [28],  The  experimental  stress  in  Figure  2.4  (a)  is 
calculated  from  the  measured  force  and  the  original  cross- 
sectional  area;  the  model  stress  is  the  average  stress  of  the 
stress  distribution  curve  in  Figure  2.4  (b)  ;  and  the  deformation 
is  the  mean  value  on  a  gage  length  of  35  mm.  The  stress 
distributions  in  Figure  2.4  (b)  were  calculated  from  the 
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Prior  to  Loading 


65  %  of  Ultimate  Load 


Figure  2.3:  Progressive  Cracking  of  Concrete  under  Uniaxial 

Compressive  Loading  [22]. 
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Figure  2.4:  (a)  Stress-Deformation  Relation  and  (b)  Calculated 

Stress  Distribution  under  Tensile  Loading  [28] . 
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measured  deformation  distributions,  the  modulus  of  elasticity 
determined  from  equilibrium,  and  the  stress-deformation 
relations  determined  from  the  average  non-elastic  effects  during 
unloading  of  the  material  adjacent  to  the  crack  faces.  The 
following  conclusions  can  be  drawn  from  Fig.  2.4:  micro-cracking 
in  tension  occurs  at  less  than  60  percentage  of  the  peak  load; 
the  stress-deformation  relation  is  linear  up  to  90  percent  of 
the  peak  load;  for  post  peak  loads  there  is  a  rapid  decrease  in 
stress  to  about  60  percent  of  that  at  the  peak  load  and 
increasing  nonlinear  deformation.  The  average  stress  at  the  peak 
load  is  about  10  percent  less  than  the  maximum  stress  at  the 
terminus  of  the  strain  softening  zone. 

By  using  moire  interferometry  on  a  dog  bone  specimen 
subject  to  tensile  stress,  Cedolin  et  al.  [26]  found  that 
concrete  fails  at  a  low  average  strain  but  that  locally  high 
strains  exist  in  the  specimen  due  to  a  non-uniform  strain 
distribution.  A  meaningful  stress-deformation  relation  is  thus 
difficult  to  obtain  from  a  direct  tension  test  since  the 
measured  average  deformation  of  the  tensile  specimen  has  little 
physical  significance  in  the  presence  of  the  localized  fracture 
zone . 


The  flexural  response  of  concrete  as  desecribed  by 
Hillerborg  et  al.  [31]  is  shown  in  Fig.  2.5. 
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Figure  2 
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.5:  Stress  Distribution  in  the  Fracture  Zone  of  a  Loaded 
Concrete  Beam  [31] . 
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2 . 3  Dynamic  Fracture  Test  Methods 


A  number  of  dynamic  test  methods  for  determining  the 
dynamic  properties  of  concrete  have  been  proposed  but  as  yet  no 
standard  test  method  has  been  established.  In  this  section,  some 
of  the  more  important  dynamic  fracture  test  methods  are  briefly 
reviewed . 

2.3.1  Drop-Weight  Test 

Nanda  and  Hannant  [32]  introduced  the  idea  of  the  "number 
of  blows  to  no  rebound"  in  order  to  reduce  the  subjectivity  and 
arbitrariness  in  the  failure  of  impacted  concrete  specimens. 
This  test  can  be  very  meaningful  when  it  is  used  in  a 
comparative  way  even  though  the  number  of  blows  depend  on  the 
drop  weight,  the  drop  height,  the  size  of  the  hammer,  its  shape 
and  the  type  of  failure.  They  found  that  plain  concrete  failed 
after  5  blows.  The  drop-weight  test  has  been  recommended  by  the 
ACI  Committee  544  to  evaluate  the  impact  resistance  of  fiber 
reinforced  concrete  [33] . 

2.3.2  Instrumented  Drop-Weight  Test 

This  test  is  a  drop-weight  test  in  which  the  striking  tup 
and  the  anvil  supports  are  instrumented  to  measure  load, 
deflection,  strain  and  energy  histories  during  the  impact  event. 
With  this  data,  static  analysis  is  used  to  compute  the  basic 
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material  properties,  such  as  ultimate  strength  and  strain  at 
peak  load,  energy  absorbed  and  fracture  toughness  at  different 
strain  rates. 

Suaris  and  Shah  [9,34]  proposed  this  method  with  a  two 
degree  of  freedom  model  for  predicting  the  amplitude  and  period 
of  inertial  loading  oscillations  for  the  tup-load  versus  time 
trace.  Mindess  and  his  co-workers  [3,15,35]  applied  this  test 
method  to  flexural  specimens  of  hardened  cement  paste,  fiber 
reinforced  concrete  and  conventionally  reinforced  concrete. 

2.3.3  Charpy  Impact  Test 

The  Charpy  impact  test  is  a  standard  impact  test 
recommended  for  metals  [36] .  The  impact  resistance  is  measured 
by  the  rebound  height  of  the  pendulum  after  impact  of  a  notched 
beam.  The  reliability  of  the  method  is,  however,  affected  by  the 
stiffness  of  the  pendulum  itself  in  measuring  the  failure  energy 
and  the  peak  load  [37].  Krenchel  [38]  and  Johnston  [39]  used 
this  method  to  measure  the  toughness  of  fiber  reinforced 
concrete . 

2.3.4  Modified  Charpy  Test 

The  modified  Charpy  test  is  a  Charpy  impact  test  in  which 
the  pendulum  and  the  anvil  supports  are  instrumented  as  in  the 
instrumented  drop-weight  test.  This  method  was  first  introduced 
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by  Hibbert  [40]  who  obtained  continuous  load  and  energy 
histories  for  plain  and  fiber  reinforced  concrete  beams 

Radomski  [41]  used  a  rotating  impact  machine  in  which  the 
impact  load  was  simulated  by  releasing  a  striker  from  a  rotating 
fly-wheel  when  it  had  attained  the  desired  velocity.  More 
recently  Shah  and  his  co-workers  [42-44]  have  conducted  several 
series  of  tests  on  plain  concrete,  mortar  and  fiber  reinforced 
concrete  specimens  using  this  method. 

2.3.5  Split  Hopkinson  Bar  Tast 

The  specimen  in  this  method  is  sandwiched  between  two 
coaxial  elastic  bars.  Tensile  pulses  are  generated  by  a  weight 
hitting  an  anvil  fixed  at  the  incident  bar.  The  force  acting  on 
the  specimen  is  calculated  by  measuring  the  pulse  at  two  points 
on  either  side  of  the  specimen  and  the  strain  in  the  specimen 
can  be  measured  by  means  of  strain  gages  mounted  on  the 
specimen . 

Bhargava  and  Rehnstrom  [45]  found  that  the  strength  of  the 
fiber  reinforced  concrete  was  about  50  percent  greater  than  the 
static  tensile  strength  for  an  observed  pulse  rise-time  of  about 
50  Usee.  Reinhardt  and  his  co-workers  [46]  developed  testing 

equipment  suitable  for  strain  rates  between  0.05  and  25/sec. 
Jawed  and  his  co-workers  [47]  generated  strain  rates  as  high  as 
about  1000/sec  and  found  that  at  a  relatively  low  strain  rate 
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(less  than  250/sec),  the  ultimate  stress  increased  with  increase 
in  strain  rate,  whereas  at  high  strain  rates  (greater  than 
250/sec),  it  reached  a  limiting  value. 

2.3.6  Dynamic  Tensila  Teat 

In  this  test,  typical  cylindrical  specimens  are  struck  at 
one  end  with  a  high  velocity  projectiles.  The  compression  wave 
generated  at  the  striking  end  is  reflected  as  a  tensile  wave 
when  it  reached  the  far  end  of  the  cylindrical  specimen  where  an 
impedance  matched  concrete  pellet  is  placed.  The  tensile  wave 
causes  the  concrete  to  spall.  The  stresses  and  strains  in 
concrete  are  found  from  the  velocity  of  the  pellet.  Birkimer  and 
Lindemann  [48]  conducted  tensile  tests  on  steel  and  nylon  fiber 
reinforced  concrete  specimens  using  this  method.  Recently,  Gran 
and  Seam  [49]  applied  this  method  to  plain  concrete  bars  and 
analyzed  the  results  using  a  smeared  crack-band  model. 

2.3.7  Explosive  Test 

This  test  consists  of  applying  an  explosive  shock  to  a 
concrete  slab.  The  spall  velocity  of  the  fragments  is  a  measure 
of  the  impact  resistance  of  the  specimen.  An  improvement  in  the 
impact  resistance  causes  a  reduction  in  the  spall  velocity  of 
the  particles  as  well  as  their  size.  This  test  is  ideally  suited 
for  situations  where  a  structure  would  be  subjected  to  rapidly 
rising  pressures  resulting  from  blasts  or  explosives. 
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2.3.8 


Constant  Strain  Rata  Test 


A  typical  test  is  carried  out  at  a  strain  rate  that  remains 
constant  throughout  the  whole  experiment.  Although  limited  in 
capacity  for  achieving  high  rates  of  loading,  conventional 
servo-controlled  machines  have  been  used  to  conduct  dynamic 
tests  at  intermediate  strain  rates.  Kobayashi  and  Cho  [50]  used 
displacement  controlled  testing  of  four-point  bend  specimens  to 
obtain  load-deflection  curves  for  polyethylene  fiber  reinforced 
concrete  beams  and  observed  that  the  loading  velocity  affected 
both  the  peak  load  carrying  capacity  and  the  corresponding 
deflection . 

2 . 4  Strain  Rate  Sensitivities  of  Mechanical  Properties 
of  Concrete 

A  number  of  attempts  have  been  made  to  determine  the 
behavior  of  concrete  under  varying  strain  rates.  Most  such 
investigations  have  shown  that  concrete  is  strain  rate 
sensitive.  Furthermore,  that  strain  rate  sensitivity  was  found 
to  differ  with  different  test  systems.  In  this  section,  the 
strain  rate  sensitivity  of  the  mechanical  properties  of  concrete 
under  compression,  tension  and  flexure  loadings  is  discussed. 

The  first  experimental  study  on  the  strain  rate  sensitivity 
of  concrete  was  conducted  by  Abrams  [51]  in  1917,  who  subjected 
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concrete  cylinders  to  impact  compression  loading.  Green  [10] 
evaluated  the  performance  of  concrete  for  varying  strain  rates 
with  the  independent  variables  being  the  cement  type,  coarse 
aggregate  type  and  shape,  curing  conditions,  grading  and  mix 
proportions,  and  the  specimen  age.  He  found  that  the  higher  rate 
sensitivities  the  greater  the  static  strength,  the  angularity  of 
aggregate  and  the  extent  of  water  curing  for  the  specimens . 

MacNeely  and  Lash  [52]  determined  the  effect  of  the  loading 
rate  on  the  flexural  strength  of  concrete  and  found  a  linear 
relationship  between  the  modulus  of  rupture,  fr/  and  the  rate  of 
loading,  G, 


fr  =  A  +  B  log  G  (2.1) 

where  A  and  B  are  constants . 

The  compression  tests  of  Atchly  and  Furr  [53]  also  showed 
that  the  compressive  strength,  absorbed  energy,  secant  modulus 
and  the  strain  at  failure  increased  with  an  increase  in  the 
strain  rate.  The  compressive  strength  and  the  energy  absorption 
reached  a  constant  limiting  value  at  high  rates  of  loading. 

Birkimer  and  Lindemann  [48]  showed  that  the  critical 
fracture  strain  energy  theory  provides  a  meaningful  fracture 
criterion  and  that  the  critical  fracture  strain  was  proportional 
to  the  strain  rate  raised  to  the  one-third  power. 
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Hughes  and  Gregory  [54]  developed  the  load  column  method  in 
which  "low  friction"  pads  was  used  to  reduce  the  platen  friction 
in  dynamic  compression  tests  on  concrete  prisms.  They  found  a 
ratio  of  dynamic  to  static  strength  of  1.9.  This  ratio  was 
independent  of  the  water-cement  ratio,  age  and  cement  content, 
but  depended  upon  the  type  of  coarse  aggregate. 

Hughes  and  Watson  [55]  found  that  the  ultimate  strains 
decreased  with  an  increase  in  the  strain  rates  for  concrete  cube 
tests  with  varying  mix  proportions  and  two  different  types  of 
coarse  aggregates  under  compressive  impact  loading  up  to  strain 
rates  of  17/sec.  This  increase  was  caused  by  the  absence  of 
creep  strains  at  high  strain  rates.  The  cause  for  the  large 
energy  requirement  in  impact  tests  was  that  more  aggregates 
failed  in  high  strain  rate  tests  whereas,  for  low  rates,  the 
crack  propagated  around  the  aggregate. 

Evans  [56]  suggested  the  velocity  of  a  growing  crack  can  be 
described  as 

V  «  Kxn  or  f'td  06  G1/(n+1)  (2.2) 

where  V  is  the  crack  velocity;  Kx  is  the  stress  intensity  factor; 
f'td  is  the  fracture  strength,  O  is  the  rate  of  stress 
application;  and  n  is  a  constant.  There  are  three  methods  of 
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20 


evaluating  the  value  of  "n"  in  Equation  (2.2).  Figure  2.6  (a) 
represents  results  from  constant  load  tests  with  the  applied 
load  causing  failure,  f'td/  plotted  against  the  time  to  failure, 
tf.  The  slope  of  that  relationship  is  1/n.  Figure  2.6  (b) 
corresponds  to  a  plot  of  the  velocity  of  the  growing  crack 
versus  Ki  and  Figure  2.6  (c)  is  the  outcome  for  a  constant  rate 
of  loading  test. 

Mindess  and  Nadeau  [57]  compared  the  slope,  n,  of  the  log  V 
versus  log  Kj  plot  with  the  slope  of  the  log  f^d  versus  log  5 

plot  for  results  from  controlled  crack  growth  tests  on  double 
torsion  specimens,  where  8  is  the  displacement  rate.  The  value 

of  "n"  was  almost  constant  in  the  case  of  mortar  but  differed  by 
a  factor  of  two  for  cement  paste,  "n"  was  also  larger  for 
compression  tests  than  for  tension  and  flexural  tests. 

Mihashi  and  Izumi  [58]  developed  a  theoretical  approach  to 
describe  the  strain  rate  results  satisfactorily  as  follows: 


where  fd  and  f3  are  the  dynamic  and  static  strengths, 
respectively;  <?d  and  are  the  dynamic  and  the  static  stress 
rates,  respectively;  and  P  is  a  material  parameter.  The 
parameter,  P,  increased  with  increased  strength  of  concrete. 


Thus  stronger  concrete  was  less  strain  rate  sensitive  than 
weaker  concrete. 

Zech  and  Wittmann  [59]  generated  strain  rates  of  about 
2/sec  by  allowing  a  missile  to  fall  on  mortar  flexural  specimens 
and  found  the  distribution  function  for  the  flexural  strength  of 
the  mortar  for  varying  strain  rates.  The  variability  of  the 
parameter,  p,  in  Equation  (2.3)  was  not  influenced  by  the  rate 

of  loading. 

Suaris  and  Shah  [9]  used  instrumented  variable  strain  rate 
tests  of  mortar  and  found  that  specimens  of  higher  flexural 
strength  were  less  strain  rate  sensitive  than  those  of  lower 
strength.  Tensile  response  was  the  most  strain  rate  sensitive, 
compressive  response  was  the  least  strain  rate  sensitive,  and 
that  flexural  response  lay  somewhere  between  those  two.  The 
constant  "n"  in  Equation  (2.2)  was  not  a  constant;  rather,  it 
decreased  with  increasing  strain  rate.  Shah  and  John  [60]  also 
point  out  that  the  log  Kj  versus  log  V  relationship  was  non¬ 
linear  especially  at  high  rates  of  loading 

Zielinski  and  Reinhardt  [61,62]  investigated  the  tensile 
stress-strain  behavior  of  mortar  and  concrete  at  high  stress 
rates  using  the  split  Hopkinson  bar  technique.  The  tensile 
strength  at  the  high  stress  rates  was  found  to  increase  markedly 
and  they  concluded  that  behavior  was  due  to  extensive  micro- 
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cracking  throughout  the  whole  volume  of  the  stressed  specimen. 
This  argument  was  supported  by  their  observation  of  higher 
ultimate  strains  at  higher  stress  rates  and  more  than  one 
fracture  surface  along  the  length  of  their  specimens.  The 
difference  between  the  impact  strength  of  concrete  and  mortar 
was  attributed  to  the  crack  arresting  mechanism  of  the 
aggregate.  For  very  rapid  loading,  cracks  were  forced  to  develop 
along  the  shorter  paths  of  higher  resistance,  through  stronger 
matrix  zones  and  also  through  some  aggregates  because  much 
energy  was  introduced  into  the  system. 

Figure  2.7  shows  the  influence  of  stress  and  strain  rates 
on  the  tensile  properties  of  concrete  such  as  tensile  strength, 
f td/ f ta »  modulus  of  elasticity,  E^/Eg,  and  strain  at  the  maximum 
stress,  eUd/£u3/  as  recommended  by  Reinhardt  [63],  The  power  a  of 

the  line  of  tensile  strength  increase  depends  on  the  concrete 
grade  as  follows: 


-  =  10  +  (2.4) 

a  z 

with  fcm  the  mean  static  cube  strength  in  N/mm2.  The  power  law 
becomes 


£td 
f  ts 


(2.5) 
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Figure  2.7:  Stress  and  Strain  Rate  Influence  on  Tensile 
Properties  of  Concrete  [63]. 
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where  Co  is  the  stress  rate  for  static  loading  with  a  value  of 
0.1  N/mm^sec.  The  power  of  the  modulus  of  elasticity  line  is 
0.016  and  that  of  the  ultimate  strain  line  is  0.020. 

Malvern  and  his  co-workers  [64]  generated  strain  rates  of 
50  to  800/sec  using  a  large  Kolsky  bar  test  system.  An  apparent 
rate  dependence  up  to  almost  twice  the  static  strength  was 
observed.  High-strength  concrete  showed  an  approximately 
logarithmic  dependence  on  the  strain  rate  at  the  maximum  stress 
while  mortar  showed  an  apparent  linear  dependence. 
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CHAPTER  3 


EXPERIMENTAL  AND  ANALYTICAL  PROCEDURES 

3 . 1  Displacement  Controlled  Dynamic  Fracture  Test  System 

A  dynamic  fracture  testing  system,  which  consisted  of  a 
high  speed  displacement-controlled  loading  system,  measurement 
system  and  data  acquisition  system,  was  constructed  for  this 
research  project.  The  displacement-controlled  dynamic  load  was 
generated  by  an  MTS  closed-loop  dynamic  testing  machine.  Figure 
3 . 1  shows  the  MTS  equipment  used  in  the  loading  system  and  the 
flow  of  commands  for  operating  the  system.  The  servo-valve 
command  is  generated  by  the  programmable  digital  wave  form 
generator  (MicroProfiler)  and  that  command  is  transmitted  to  the 
displacement  controller  (AC  Controller)  through  the  system 
controller  (MicroConsole) .  The  displacement  controller  generates 
the  inner  loop  (IL)  command  to  the  valve  controller.  The  valve 
controller  and  servo-valve  linear  variable  differential 
transformer  (LVDT)  form  an  inner  loop  to  produce  an  error  signal 
which  represents  the  difference  between  the  desired  and  actual 
transducer  output.  The  error  signal  is  conditioned  to  produce  a 
valve  command  which  controls  the  mechanical  action  of  the 
actuator.  The  programmed  command  of  the  function  generator  is 
also  used  as  the  triggering  signal  for  the  data  acquisition 
system.  The  340  1/min  (90  gpm)  servo-valve  and  380  1/min  (100 
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Figure  3.1:  MTS  Dynamic  Loading  System 


gpm)  hydraulic  service  manifold  used  in  the  system  can  generate 
a  maximum  piston  velocity  of  3.38  m/s  (133  in/sec).  An  48.9  kN 
(11  kip)  actuator  is  mounted  on  the  stroke  crosshead  and  the 
stiffness  of  the  loading  system  is  provided  by  a  245  kN  (55  kip) 
tee-slotted  table  and  two  columns. 

The  gages  and  the  devices  for  the  measurement  system  were 
carefully  selected  for  high  speed  measurement  of  70  KHz  resonant 
frequency.  The  system,  as  configured  for  testing  of  a  three 
point  notched  beam  specimen,  is  shown  in  Figure  3.2.  The  system 
consists  of:  a  dynamic  load  cell  and  a  power  unit;  three  strain 
gages  and  a  ten-channel  strain  gage  conditioner/amplifier;  and 
two  inductive  displacement  gages  and  two  oscillator/ 
demodulators . 

The  linear  capacity  of  the  dynamic  load  cell  is  22.3  kN 
(5  kips)  and  the  load  cell  is  designed  for  fast  transient  and 
repetitive  applications.  It  has  an  inertial  mass  of  14.2  g  (0.5 
oz)  to  minimize  its  inertial  effect.  The  load  cell  consists  of 
two  thin  quartz  discs  loaded  in  compression  and  sandwiched 
between  hardened  steel  cylindrical  members . 

The  foil  strain  gages,  5  x  1.5  mm,  are  located  along  the 
crack  path  to  monitor  the  transient  strains  associated  with  a 
propagating  crack.  Inductive  displacement  gages  or  non-contact 
capacitance  gages  of  4.1  mm  (0.16  in)  maximum  capacity  are  used 
for  transient  crack  opening  displacement  measurement. 
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Figure  3.2:  Dynamic  Measurement  System  and  a  Three-Point  Bend 
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For  the  data  acquisition  system,  three  digital-storage 
oscilloscopes  and  a  personal  computer  were  used.  The  data  in  the 
oscilloscopes  was  subsequently  transferred  to  the  personal 
computer  for  processing,  interpretation  and  permanent  storage. 
Each  oscilloscope  has  two  channels  with  each  channel  capable  of 
storing  up  to  1,024  data  per  sampling  within  a  minimum  of  0.001 
sec  and  a  maximum  of  2  sec. 

The  configuration  of  the  system  for  testing  of  CLWL-DCB 
specimens  is  described  in  Chapter  6.  Details  of  the  dynamic 
loading  system  and  the  instrumentation  for  those  tests  were  the 
same  as  described  for  the  three  point  bend  specimens. 

3 . 2  Impact  Loading  Test  System 

The  strain  rates  achieved  in  the  displacement  controlled 
three  point  beam  tests  did  not  reach  the  maximum  values  desired 
for  this  research.  Further,  the  numerical  analyses  indicated 
relative  energy  contribution  (Fig.  4.10)  significantly  different 
from  this  typically  associated  with  impact  loading  results  (Fig. 
5.10).  Accordingly,  an  impact  loading  test  set-up  was  also 
developed  in  an  environment  where  moire  interferometry  could 
also  be  used  to  measure  crack  opening  displacements. 

A  special  testing  machine  was  constructed  in  the  optics 
laboratory  of  the  Department  of  Mechanical  Engineering  for  the 
impact  loading  and  moire  interferometry  studies.  The  impactor 
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weighed  98.8  N  (22.2  lb)  and  the  maximum  height  of  drop  was  1.0 
m  (39.4  in)  .  The  impactor  was  dropped  from  a  height  of  0.4  m 
(15.7  in).  In  the  same  manner  as  in  the  displacement  controlled 
tests  the  concrete  specimens  were  instrumented  with  a  load  cell 
at  the  contact  point  and  three  strain  gages  located  at  the  same 
points  as  shown  in  Fig.  3.2.  Figure  3.2(b)  is  the  test  set-up 
for  the  impact  loading  three  beam  tests. 

The  basic  principles  of  Moire  interferometry  were  described 
in  Chapter  Four,  Part  I  of  this  contract  report.  A  specimen 
grating  of  600  lines/mm  and  a  reference  grating  of  1200  lines/mm 
were  used  for  moire  interferometry.  The  specimen  grating  was 
illuminated  by  an  argon  ion  laser  at  an  oblique  incident  angle 
of  about  17°  where  the  first  order  diffracted  light  provided  a 
fringe  multiplication  of  the  specimen  grating  and  interfered 
with  the  reference  grating.  For  dynamic  fracture  analysis,  the 
transient  Moire  interferometry  fringes  were  recorded  by  an 
IMACON  790  high  speed  camera  and  1/5T  plug-in  with  a  framing 
rate  of  100,000  frames  per  second  and  an  exposure  time  of  2 
microseconds.  This  camera  was  specially  configured  to  provide 
four  large  frames  (323  x  18  mm)  .  The  larger  frames  were 
necessary  to  resolve  the  dense  Moire  fringes.  A  high  speed 
electric  shutter  with  an  exposure  time  of  7  milliseconds  was 
inserted  into  the  optical  path  to  shield  the  image  tube  of  the 
camera  from  prolonged  exposure  to  the  argon  ion  laser  light 
source . 
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3 . 3  Experimental  Test  Program 


The  effects  of  geometry  and  overall  stress  state  have  been 
investigated  by  conducting  multiple  tests  on  two  different  type 
tests  of  fracture  specimens:  three-point  bend  specimens  and 
CLWL-DCB  specimens.  The  same  cement,  sand,  aggregate  and 
aggregate  gradation  were  used  for  all  three  point  bend  and  CLWL- 
DCB  specimens  so  that  interpretations  of  the  significance  of  the 
results  could  be  concentrated  on  strain  rate  and  geometrical 
sensitivity  issues  rather  than  material  variation  issues. 
Specimens  were  prepared  using  the  same  procedures  as  those 
developed  for  prior  static  tests  [13] .  Tables  3.1  and  3.2  show 
the  average  gradations  for  coarse  and  fine  aggregates, 
respectively.  The  cement  was  a  type  III,  high  early  strength 
Portland  cement,  and  the  sand  and  the  gravel  were  obtained 
locally.  The  same  mix  proportions  were  used  for  all  specimens 
and  those  proportions  are  shown  in  Table  3.3. 
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Table  3.1:  Average  Gradation  of  Coarse  Aggregate 


0.75 

tr 

0.37 

5" 

0.25 

// 

No. 

4 

No  . 

8 

1  % 

0 

0 

0 

0 

72 

100 

Table  3.2:  Average  Gradation  of  Fine  Aggregate 


Sieve 

Size 

No .  4 

No.  8 

No.  30 

No.  50 

No. 

100 

% 

1.50 

15.79 

28.87 

43.88 

78.54 

97.50 

Table  3.3:  Concrete  Mix  Proportion  by  Weight 


Cement 

Sand 

Gravel 

Water 

1.0 

2.5 

2.0 

0.4- 

0.5 

The  specimens  were  cast  in  a  horizontal  position  in  oil 
lubricated  steel  forms  and  steel  cylinders.  For  the  CLWL-DCB 
specimens,  those  forms  contained  inserts  screwed  to  them  to 
create  holes  for  the  split  pins  and  the  starter  cracks.  Each 
specimen  was  compacted  with  a  mechanical  form  vibrator  that 
moves  across  the  top  of  the  specimen.  Particular  care  was  taken 
to  avoid  segregation  of  the  ingredient  materials.  Control 
cylinders  were  cast  with  each  batch  of  fracture  specimens  and 
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each  control  cylinder  fabricated  in  three  layers  and  each  layer 
compacted  with  a  one-inch  diameter  electric  immersion  vibrator. 

All  specimens  and  their  companion  cylinders  were  cured  in 
the  same  manner.  Immediately  after  casting,  the  specimens  and 
cylinders  were  covered  with  polyethylene  sheets  and  stored  at 
room  temperature.  The  forms  were  stripped  at  the  age  of  two 
days;  steel  inserts  carefully  removed  from  the  specimens;  and 
then  both  specimens  and  cylinders  immediately  placed  in  a  moist 
room.  Each  specimen  was  removed  from  the  moist  room  two  days 
before  testing.  After  it  had  dried  for  one  day,  the  surface  of 
the  specimen  was  polished  and  all  necessary  gages  bonded  to  that 
surface . 

3 . 4  Numerical  Analysis 

3.4.1  Dynamic  Finite  Element  Analysis 

The  finite  element  method  used  in  the  research  was  derived 
from  the  regular  four-node  quadrilateral  element  using  the 
Hermitian  interpolation  function  [65]  and  two-dimensional 
element  formulations. 

For  dynamic  analysis,  Newmark's  (3  method  is  applied  to 
solve  the  equation  of  motion, 

[K]  {u}  +  [M]  {ii}  =  {P},  (3.1) 
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approximated  by  the  finite  element  model,  where  [K]  and  [M]  are 
the  stiffness  and  mass  matrices,  respectively;  {u}  and  {u}  are 
the  displacement  and  acceleration  vectors,  respectively,  of  the 
finite  element  assemblage;  and  {P}  is  the  load  vector. 

By  using  Newmark's  (3  method  as  a  direct  integration  method, 

all  solution  techniques  employed  in  static  analysis  can  be  used. 
Another  advantage  of  using  Newmark's  {3  method  as  an  implicit 

scheme  is  its  unconditional  stability  since  the  implicit  scheme 
computes  the  displacement  vector  directly  from  the  equation  of 
motion  at  the  current  time  step  while  the  explicit  scheme 
evaluates  the  current  displacement  vector  using  the  acceleration 
vector  computed  at  the  previous  time  step. 

3.4.2  Application  of  Crack  Closure  Stress 

To  simulate  the  softened  stiffness  of  the  fracture  process 
zone,  crack  closure  forces  were  applied  at  released  nodes.  The 

equivalent  nodal  force  for  the  crack  closure  stress 

corresponding  to  a  given  crack  opening  displacement,  to,  the 

crack  closure  force  vector,  Fccs,  of  a  released  crack  tip  element 
i  was  calculated  as  follow: 

Fccs  =  {tO} 

=  J  [B]T[C„]  {tO}  dV 

Vl 
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(3.2) 


=  |  [B] T { fcca }  dV 

Vl 

=  t  J  [B ] T { fCcs )  dA 
*1 


where  [k^]  is  the  element  stiffness  matrix  for  the  crack  opening 
displacement;  Vi,  Ai  and  t  are  the  volume,  area  and  thickness  of 
the  released  element  i;  and  [B]  is  the  strain-displacement 
matrix  of  the  released  element.  The  nodal  stress  vector,  (fees)/ 
is  calculated  from  the  assumed  constitutive  relation,  [C<n]  , 

between  the  crack  closure  stress  and  crack  opening  displacement. 
The  direction  of  the  crack  closure  force  vector  is  opposite  to 
that  of  the  crack  opening  displacement. 

3.4.3  Modelling  of  Continuous  Crack  Propagation 

Crack  propagation  is  simulated  by  releasing  constrained 
nodes.  For  a  large  crack  extension  in  a  short  time  step,  this 
procedure  releases  a  large  unbalanced  potential  energy. 

This  discrete  change  is  moderated  by  a  gradual  movement  of 
the  physical  crack  tip  in  element,  i,  as  shown  in  Figure  3.3. 
Here  the  distance  of  the  physical  crack  tip  is  assumed  to  be 
linearly  dependent  on  the  stress  at  the  mathematical  crack  tip 
(node  j+1)  as: 


(3.3) 
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Figure  3.3:  Numerical  Simulation  of  Continuous  Crack  Propagation 
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where  |3  is  the  damage  rate  for  the  crack  tip  element  i;  s  is  the 
distance  of  the  physical  crack  tip  from  the  released  node  j;  Li 
is  the  mesh  size  of  the  crack  tip  element  in  the  crack 
propagation  direction;  ft  is  the  tensile  strength;  and  f0  and  fj+i 
are  the  stresses  at  the  mathematical  crack  tip  node  j+1  when  the 
node  j  is  released  and  at  the  current  time,  respectively. 

To  simulate  effects  for  the  uncracked  portion  of  the  crack 
tip  element  i,  a  crack  tip  force,  Ftip,  is  also  applied  at  the 

released  node  j  of  the  crack  tip  element.  That  crack  tip  force 
is  determined  from  the  damage  rate  [5  and  the  response  force  for 

the  crack  tip  node  j  when  that  node  is  released  as  follows: 

Ftip  =  (3F°ip  (3.4) 

where  Ftip  is  the  crack  tip  force  required  to  simulate  the 
stiffness  of  the  uncracked  portion  of  the  crack  tip  element  i 
and  F°ip  is  the  response  force  of  node  j  when  that  node  is 
released. 

For  the  cracked  portion  of  the  crack  tip  element  i,  crack 
closure  force  effects  are  calculated  in  a  similar  manner  as 
follows : 

Fcci  -  <H»  F=C,  <3-5> 
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where  F*'^  is  the  crack  closure  force  for  the  cracked  portion  of 
the  crack  tip  element  i  acting  on  the  released  crack  tip  node  j 
and  F^cs  is  the  crack  closure  force  for  the  crack  tip  element  i 
calculated  from  the  crack  opening  displacement  at  node  j  as 
described  in  Section  3.4.2.  Therefore  the  total  force  applied  to 
the  released  node  j  of  the  crack  tip  element  i  becomes 


F  tip  + 


D 

CCS 


+  F 


i-1,  j 

CCS 


(3.6) 


where  F^.C3' 3  is  the  fraction  of  the  total  crack  closure  force 
associated  with  the  released  element  i-1  and  acting  on  the 
node  j  as  the  result  of  the  crack  opening  displacement  of  the 
node  j . 

3 . 4  Determination  of  Mechanical  Properties  of  Concrete 

The  mechanical  properties  of  concrete  under  dynamic  load 
were  investigated  using  a  hybrid  experimental-numerical 
technique.  First  the  elastic  properties  of  the  test  specimens 
were  determined  from  the  experimental  data  using  the  measured 
responses  prior  to  micro-crack  propagation.  Then  the  fracture 
properties  were  optimized.  Figure  3.4  shows  the  flow  chart  of 
the  optimization  procedure  for  determining  the  elastic  and 
fracture  properties. 
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Figure  3.4:  Flow  Chart  for  Determination  of  Mechanical 

Properties 
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3.4.1  Hybrid  Exp«ria«ntal-Num*rical  Technique 

The  hybrid  experimental-numerical  technique  [66,67] 
consists  of  two  phase  analyses  using  "generation"  and 
"application"  modes.  In  the  generation  mode,  a  characteristic 
form  is  assumed  for  the  crack  closure  stress  function  and  the 
parameters  quantifying  that  form  are  determined  through  a 
numerical  technique  using  boundary  conditions  observed 
experimentally.  In  the  application  mode,  the  numerical  technique 
is  used  to  generate  results  for  an  extending  crack.  Those 
results  are  compared  with  the  experimental  results.  Postulated 
boundary  conditions  and  constitutive  equations  are  then  modified 
and  the  procedure  repeated  until  the  numerical  and  experimental 
results  are  in  agreement  with  each  other. 

In  this  research,  the  load  history  was  used  as  an  input 
boundary  condition  for  a  two  dimensional  dynamic  finite  element 
analysis  and  strain  histories  were  generated  numerically.  Those 
numerically  generated  strain  histories  were  then  compared  with 
the  strain  histories  measured  experimentally  to  optimize  the 
elastic  and  fracture  constitutive  relationships. 

This  hybrid  experimental-numerical  method  provides 
reasonable  results  since  it  does  not  require  unrealistic 
assumptions,  such  as  a  constant  stress  state  for  a  direct  tensile 
test  or  beam  theory  for  a  modulus-of-rupture  test.  Also,  the 
method  can  accomodate,  in  theory,  variable  mechanical  properties 
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throughout  the  fracture  specimen  as  predicted  from  available 
data  on  the  strain  rate  sensitivities  of  tension,  compression 
and  flexure  tests  [60,68,69]. 

3.4.2  Elastic  Propsrtias 

With  the  hybrid  experimental-numerical  technique,  the 
mechanical  properties  required  for  the  concrete  can  be 
determined  by  numerical  analysis.  The  elastic  numerical  analysis 
requires  the  following  parameters:  moduli  of  elasticity  in 
tension  and  compression;  and  tensile  strength  or  maximum  tensile 
strain . 

Moduli  of  elasticity  were  determined  by  matching  the 
initial  part  of  the  numerical  and  experimental  strain  histories. 
Since  analysis  based  on  strain  gage  data  is  very  sensitive  to 
the  modulus  of  elasticity,  it  is  necessary  to  use  different 
moduli  of  elasticity  for  tensile  and  compressive  elements,  where 
the  tensile  and  compressive  elements  are  defined  from  the  rate 
of  volume  change.  Once  the  moduli  of  elasticity  are  determined, 
the  crack-initiation  time  can  be  determined  as  the  time  when  the 
numerical  strain  no  longer  matches  the  experimental  strain.  The 
tensile  strength  was  calculated  numerically  from  the  maximum 
tensile  stress  at  the  crack-starting  time.  Although  this  tensile 
strength  is  dependent  on  the  mesh  size,  that  dependence  can  be 
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Figure  3.5:  FPZ  Constitutive  Relation 
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overcome  by  taking  the  maximum  aggregate  size  as  the  mesh  size 
because  the  experimental  value  for  the  strength  of  concrete  also 
depends  on  the  maximum  aggregate  size. 


3.4.3  Fracture  Properties 

The  fracture  behavior  of  the  cohesive  crack  model  is 
specified  using  a  softening  curve.  The  softening  curve  is 
described  using  the  following  parameters:  the  tensile  strength; 
the  shape  factor;  and  the  critical  crack  opening  displacement, 
as  shown  in  Figure  3.5. 


In  this  research,  the  softening  curve  was  assumed  as 


f  CC3 


(3.7) 


where  fees  is  the  crack  closure  stress,  which  represents  the 
remaining  stiffness  of  the  fracture  process  zone;  ftd  is  the 
dynamic  tensile  strength;  co  and  (0C  are  the  crack  opening 

displacement  (COD)  and  the  critical  crack  opening  displacement 
(CODc) ,  respectively;  and  a  is  the  shape  factor.  Therefore,  the 
fracture  properties  are  specified  with  the  shape  factor,  a,  and 
the  CODc  values.  Those  two  parameters  are  optimized  by  comparing 
the  strain  histories,  which  are  generated  numerically  from  the 
experimental  load  history,  with  the  strain  histories  observed 
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experimentally  after  the  crack-starting  time.  The  fracture 
energy  density  Gp  of  the  assumed  softening  curve  is  calculated  as 

Wo 

Gp  =  J  fees  d(0  =  ftd  0>c  (3.8) 

o  l  +  a 

It  is  anticipated  that  the  fracture  energy  density  Gp  is  a 
material  property  independent  of  the  strain  rate  and  the 
specimen  geometry.  For  static  loading,  the  fracture  energy 
density,  Gp,  and  the  stiffness  in  the  FPZ  were  chosen  to  match 
those  provided  by  the  corresponding  three  segment  model  of  [68]. 
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CHAPTER  4 


DISPLACEMENT  CONTROLLED  THREE-POINT  BEND  TESTS 

4 . 1  Experimental  Set-Up 

Three-point  bend  tests  were  used  in  the  first  phase  of  this 
research  since  such  tests  are  relatively  simple  to  perform  and 
the  measured  data  can  be  easily  interpreted.  Figure  4.1  shows 
the  dimensions  of  the  three-point  bend  specimen  used,  and  the 
arrangement  of  gages  and  fixtures.  An  initial  crack  of  a  12.7  mm 
(0.5  in)  length  and  3.2  mm  (0.125  in)  width  was  precast  into  the 
specimen  using  a  removable  steel  plate.  That  starter  crack 
insured  straight  crack  growth. 

An  actuator  mounted  on  the  crosshead  was  used  to  generate 
external  load  by  vertical  motion  of  its  piston  as  described  in 

Section  3.1.  A  pre-load  of  about  22.2  N  (5  pounds)  was  applied 

\ 

for  the  stable  positioning  of  the  specimen,  load  block  and 
supports.  The  piston  was  then  accelerated  to  its  programmed 
velocity  and  moved  at  a  predetermined  constant  velocity 
thereafter  until  the  visible  crack  tip  reached  more  than  90 
percent  of  the  beam  depth. 

The  signal  from  the  function  generator  that  the  test  was 
starting  triggered  the  three  oscilloscopes  allowing  the  output 
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Figure  4.1:  Single-Edge  Notched  Three-Point  Bend  Specimen 


voltages  from  the  six  gages  to  be  recorded  in  the  storage  of  the 
oscilloscopes.  The  response  load,  including  the  inertial  effect 
of  the  steel  load  block,  was  measured  by  a  dynamic  load  cell 
fixed  on  the  top  of  the  load  block.  Two  non-contact  capacitance 
gages,  which  were  mounted  on  either  side  of  the  load  block, 
measure  the  beam  deflection  or  load-line  displacement  relative 
to  MTS  machine  table.  Three  foil  strain  gages,  5  x  1.5  mm,  were 
bonded  at  12.7  mm  (0.5  in)  offsets  from  the  center  line  of  the 
beam.  Despite  that  offset,  those  gages  provided  information  not 
only  on  the  strain  histories  in  the  fracture  process  zone  (FPZ) 
but  also  the  time  when  the  micro-crack  tip  of  the  FPZ  passed 
those  gages.  That  was  possible  because  the  speed  at  which  the 
unloading  compression  stress  waves  from  a  moving  crack  tip 
propagated  laterally  towards  the  gage,  was  very  fast  compared  to 
the  relatively  slow  velocity  of  the  propagating  crack.  That 
crack  velocity  was  less  than  five  percent  of  the  elastic 
longitudinal  wave  speed. 

From  the  bend  test,  five  experimental  results  were 
obtained:  the  load  history;  the  load-line  displacement  history 
(the  average  reading  of  two  inductive  gages)  ;  and  the  strain 
histories  for  three  different  depths  at  which  gages  were 
attached.  Of  those  results,  the  load-line  displacement  history 
was  the  least  reliable  result  because  values  include  effects 
associated  with  penetration  of  the  load  and  the  supports  into 
the  speciment  during  the  test. 
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4 . 2 


Numerical  Analysis 


Figure  4.2  shows  a  typical  finite  element  model  used  in 
Lng  the  data  lor  a  three-point  bend  specimen.  The 

_ ...  _..ical  properties  of  the  specimens  were  optimized  by 

comparing  the  strain  histories  measured  experimentally  with 
those  generated  numerically  from  the  measured  load  history  and 
the  assumed  mechanical  properties.  The  validity  of  the  optimized 
mechanical  properties  was  then  verified  by  comparing  the 
measured  load-line  displacement  history  with  the  numerically 
generated  one . 

Initiation  of  dynamic  crack  propagation  at  the  precast 
notch  tip  was  complicated  by  pre-existing  micro-cracks  generated 
bv  shrinkage  stresses.  Through  a  trial-and-error  process,  this 
ect  was  found  to  be  best  modelled  by  reducing  the  stiffness 
f  the  crack  tip  element  and  delaying  the  initiation  time  of 
crack  extension  as  that  given  by  the  strain  history  of  gage 
SG01 . 


The  strain  rate  sensitivity  of  the  dynamic  tensile 

»  ,  ,  , 

strength,  ffcd,  was  assumed  equal  to  the  strain  rate  sensitivity 

of  the  maximum  impact  load.  The  maximum  static  load  used  to 
normalize  that  strain  rate  sensitivity  was  the  maximum  load  for 
the  Group  5,  tests,  which,  as  apparent  from  Table  4.1,  were  the 
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Figure  4.2:  Finite  Element  Model  of  Single-Edged  Notched  Three-Point  Bend  Specimen 


slowest  loading  tests.  The  fracture  properties,  a  and  CODc,  as 
well  as  the  tensile  and  compressive  moduli  of  elasticity,  Etd  and 
Ecd/  respectively,  were  obtained,  as  described  in  Chapter  3,  via 
an  inverse  process  in  which  these  properties  were  optimized  to 
fit  the  measured  strain  histories. 

As  in  all  inverse  procedures,  the  uniqueness  of  the 
constitutive  relation  for  the  FPZ  and  the  tensile  and 
compressive  moduli  of  elasticity  are  not  guaranteed.  The 
computed  results  were  assumed  correct  within  the  region  where 
the  computed  strains  at  the  three  strain  gage  locations  and  the 
computed  load-line  displacement  agreed  with  their  measured 
counterparts.  Another  validation  check  involved  computing  the 
kinetic,  strain  and  fracture  energies  separately  and  then 
comparing  the  sum  of  these  three  quantities  with  the  total  input 
work  at  each  instance  of  crack  extension.  The  final  validation 
check  involved  the  internal  consistency  of  the  computed  results 
and  their  convergence  to  their  known  static  counterparts  at  low 
strain  rates.  In  spite  of  the  trial-and-error  process  used  in 
this  analysis,  these  strong  correlations  between  the  several 
different  computed  and  measured  data  left  little  reason  to 
suspect  that  the  computed  material  properties  deviated  from  the 
corresponding  measured  results. 
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4 . 3 


Results 


4.3.1  Experimental  Results 

Five  groups  of  specimens,  with  five  identical  specimens  in 
each  group,  were  loaded  at  five  different  programmed  actuator 
displacement  rates  which  yielded  failure  times  of  2,  5,  25,  100, 
and  1000  msec.  Of  the  five  tests,  two  test  results  which  best 
represented  the  responses  of  each  of  those  five  groups,  were 
used  in  subsequent  numerical  analyses.  The  test  conditions  for 
each  of  these  two  specimens  are  listed  in  Table  4.1. 


Table  4 . 1 

Three-Point  Bend  Test  Data 


Group 

No. 

Test 

No. 

Pmax 

<kN) 

Time  to 
Pmax 
(msec) 

SGOl 

max. 

Strain  (p£) 
SG02 

max. 

SG03 

min. 

Time  to  Failure 
(msec)  tF 

1 

11 

7.34 

1.5 

43 

115 

-80 

2.0 

12 

7.43 

1.4 

68 

144 

-93 

2.0 

2 

21 

5.56 

3.8 

47 

144 

-84 

4.6 

22 

5.29 

3.5 

31 

123 

-64 

4.3 

3 

31 

4.25 

18.6 

40 

123 

-78 

21.9 

32 

4.29 

18.6 

46 

102 

-71 

21.9 

4 

41 

3.88 

74.8 

36 

84 

-62 

84.0 

42 

3.92 

71.7 

49 

101 

-67 
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Typical  load  and  strain  histories  for  a  specimen,  Test  11, 
loaded  to  failure  in  2  msec  are  shown  in  Figure  4.3.  The  time  at 
which  the  FPZ  tip  was  taken  as  passing  a  given  gage  position,  is 
indicated  by  a  symbol  on  the  response  for  that  gage.  The  strain 
peaks  for  strain  gages  SG02  and  SG03  in  Figure  4.1,  were  taken 
as  the  time  at  which  the  FPZ  tip  passed  those  gages.  However, 
such  was  not  the  case  for  the  first  strain  gage,  i.e.  SG01  in 
Figure  4.1,  since,  as  discussed  previously,  numerical 
experimentation  demonstrated  that  locked-in  shrinkage  effects 
altered  the  behavior  at  that  location.  The  first  plateau  in  the 
SGOl  record  indicated  that  dynamic  microcrack  tip  extension 

I 

initiated  at  1.10  msec  and  at  a  stress  lower  than  The 
continued  slight  rise  and  the  subsequent  plateau  in  the  SGOl 
strain  history  represent  continuous  increase  in  strain  in  the 
crack  tip  element  despite  the  onset  of  micro-crack  propagation 
at  1.10  msec.  In  contrast,  the  passage  of  the  micro-crack  tip 
past  SG02  and  SG03  is  characterized  by  a  precipitous  drop  in 
strains.  Average  crack  velocities  were  64  m/sec  (2,500  in/sec) 
between  strain  peaks  for  SGOl  and  SG02  and  154  m/sec 
(6,060  in/sec)  between  strain  peaks  for  SG02  and  SG03 . 

Figure  4.4  shows  the  load  and  the  three  strain  histories 
for  the  specimen,  Test  52,  loaded  to  failure  in  1000  msec.  The 
average  crack  velocities  were  0.060  m/sec  (2.38  in/sec)  between 
the  strain  peaks  for  SGOl  and  SG02  and  0.233  m/ssec  (9.17 
in/sec)  between  the  strain  peaks  for  SG02  and  SG03 . 
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Figure  4.3:  Load  and  Strain  Histories  of  Test  11 
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Figure  4.4:  Load  and  Strain  Histories  of  Test  52 


Figure  4.5  shows  the  load  versus  load-line  displacement 
histories  for  the  five  different  strain  rate  loadings.  The 
increase  in  strength,  and  to  a  lesser  extent  initial  stiffness, 
with  increasing  strain  rate  is  obvious. 

4.3.2  Numerical  Results 

The  experimental  results  for  two  tests  from  each  of  the 
five  groups  of  tests  were  analyzed  numerically  using  the 
measured  loading  history  and  the  postulated  constitutive 
relation  for  the  FPZ,  i.e.  Equation  (3.8),  as  input  boundary 
conditions.  This  propagation  analysis  was  preceded  by  an 
extensive  sensitivity  analysis  of  the  variables  involved  in  this 
inverse  process.  The  particular  form  of  the  continuous  function 
for  the  constitutive  equation  of  the  FPZ  was  based  on  the 
previously  obtained  three  segment  model  of  the  FPZ  [68,69].  Use 
of  that  static  model  greatly  reduced  the  amount  of  trial-and- 
error  analysis  involved  in  matching  the  computed  and  measured 
crack  extension  and  the  three  strain  gage  histories. 

Figure  4 . 6  shows  the  computed  crack  extension  histories  for 
the  specimens  which  were  loaded  to  failure  in  4.6  and  84  msec, 
respectively.  The  slow  crack  velocities  during  the  initial  phase 
of  dynamic  crack  propagation  represent  the  influence  of  the  pre¬ 
existing  damage  region,  at  the  precast  notch  tip,  caused  by 
shrinkage  effects.  The  symbols  on  Fig.  4.6  indicated  measured, 
as  opposed  to  computed  results . 
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Figure  4.5:  Load  versus  Load-Line  Displacement 
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SG01 


Figure  4.6:  Crack  Extension  Histories  of  Tes 


Figures  4.7  and  4.8  show  the  measured  and  computed  strain 
histories  at  the  three  strain  gage  locations  for  Group  2  and  4 
tests.  Excellent  agreements  between  the  computed  and  measured 
strains  are  noted  prior  to  the  arrival  of  the  microcrack  tip  of 
the  FPZ  at  the  location  of  SG02 .  The  subsequent  decrease  in  the 
measured  strains  is  also  modelled  relatively  well  by  the 
numerical  simulation. 

Figure  4 . 9  shows  the  measured  and  computed  load  versus 
displacement  relationships  for  Tests  21  and  41.  Excellen. 
agreements  between  the  measured  and  numerical  results  are  noted. 
Much  of  the  strain  softening  portion  of  the  loading  history  was 
not  simulated  in  Figure  4.9  due  to  limits  on  the  amount  of 
available  computer  time. 

Figure  4.10  shows  the  variations  in  energy  partitions  with 
crack  propagation  for  Tests  21  and  41.  The  sums  of  the  three 
energies,  i.e.  strain,  kinetic  and  fracture  energies,  were  at 
the  worst  within  3  percent  of  the  corresponding  input  work 
during  the  dynamic  fracture  process.  It  is  noteworthy  that, 
unlike  similar  impacted  metallic  three-point  bend  specimens,  the 
kinetic  energy  is  negligible.  Because  of  the  lack  of  substantial 
kinetic  energy,  the  two  broken  halves  of  the  bend  specimens 
remained  attached  under  this  displacement  controlled  impact 
loading  condition.  The  good  to  excellent  agreement  between  the 
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Experimental  Data 


Figure  4.7:  Experimental  and  Computed  Strain  Histories  of  Test  21 


(N>|)  pool 


Computed  Strain  Histories  of  Test  41 
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measured  and  computed  values  in  the  foregoing  comparison 
provides  additional  verification  of  the  numerical  analysis.  The 
final  material  properties  gleamed  from  the  entire  numerical 
analysis,  two  per  group  for  the  five  groups,  are  summarized  in 
the  following  section. 

4.3.3  Strain  Rata  Sensitivity  of  Mechanical 
Properties 

The  computed  strain  rate  sensitivities  of  the  tensile 
strength  and  the  tensile  and  compression  moduli  of  elasticity 
are  shown  graphically  in  Figures  4.11  and  4.12,  respectively. 
Symbols  indicate  computed  results  and  the  bilinear  curves 
represents  the  best  fit  to  those  results.  The  higher  dynamic 
tensile  strength,  ftd,  at  higher  strain  rates  is  consistent  with 
the  observed  increase  in  trans-aggregate  failures  on  the 
fracture  surfaces.  In  contrast,  the  fracture  surfaces  of 
statically  loaded  specimens  exhibited  interface  failures  around 
the  aggregates  due  to  the  lower  interface  strength  between  the 
aggregates  and  the  cement  matrix. 

The  two  pseudo-static  results  control  the  horizontal 
portions  of  the  bilinear  curves  and  from  Figure  4.12  Et  is 
predicted  to  be  0.83  Ec.  The  results  for  two  representative  tests 
from  each  of  the  five  groups  and  the  associated  numerical 
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Figure  4.11:  Strain  Rate  Sensitivity  of  Tensile  Strength 
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Figure  4.12:  Strain  Rate  Sensitivity  of  Moduli  of  Elasticity 


analyses  are  summarized  in  Table  4.2.  Average  strains  in  the 
initial  crack  tip  element  prior  to  onset  of  micro-crack  tip 
propagation  were  used  to  compute  the  strain  rates  as  shown  in 
column  3.  The  damaged  region  which  was  induced  by  shrinkage 
stress  at  the  precast  notch  tip  during  concrete  curing  was 
modeled  by  using  the  reduced  modulus  of  elasticity,  Etsm  shown  in 
column  6. 

Figure  4.13  shows  the  variations  in  the  constitutive 
equation,  which  relates  the  crack  closure  stress  and  the  COD  in 
the  FPZ ,  with  strain  rate.  The  crack  closure  stress  shape 
factor,  a,  (Fig.  3.5  and  Eq.  3.7)  is  related  to  the  slope  of  the 
constitutive  relation  in  Figures  4.2  and  4.13.  Values  for  a  are 
listed  in  column  (8)  of  Table  4.2.  Also  shown  for  comparison  is 
the  three-segment  model  which  was  obtained  previously  for  static 
crack  extension  [68]. 

Figure  4.14  shows  that  the  fracture  energy  density,  Gp,  is 
strain  rate  insensitive.  Its  value  remained  constant  in  spite  of 
the  application  of  loading  rates  that  markedly  changed  the 
tensile  strength  and  modulus  of  elasticity  of  the  concrete. 

4  .  4  Discussion 

The  increasing  tensile  strength  with  increasing  loading 
rate  indicates  a  reduction  in  the  microcracked  region,  which 
precedes  the  propagating  crack  tip,  with  increasing  strain  rate. 
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Table  4.2 

Strain  Rate  Sensitivity  of  Mechanical  Properties 


Group 

No. 

Specimen 

No. 

^(sec-1) 

Ecd 

(GPa) 

E* 

(GPa) 

Elm 

(GPa) 

* 

(\$a) 

a 

Gf 

(N/m) 

1 

Oil 

0.231 

58.6 

55.2 

12.4 

7.86 

1.6 

154 

012 

0.252 

58.6 

55.2 

19.3 

8.00 

1.6 

156 

2 

021 

0.108 

49.6 

42.7 

6.2 

5.98 

1.8 

154 

022 

0.0877 

51.7 

43.4 

9.0 

5.67 

1.8 

154 

3 

031 

0.0189 

43.4 

41.4 

8.3 

4.56 

2.0 

155 

032 

0.0158 

44.8 

44.8 

13.8 

4.61 

2.0 

154 

4 

041 

0.0046 

44.8 

40.0 

9.0 

4.17 

2.2 

149 

042 

0.00373 

44.8 

37.9 

15.2 

4.21 

2.2 

151 

5 

051 

0.00053 

41.4 

34.5 

6.9 

3.66 

2.5 

144 

052 

0.00045 

41.4 

34.5 

13.8 

3.79 

2.5 

144 

£:  average  strain  rate  computed  at  notch  tip  element  prior 

to  onset  of  microcrack  propagation. 

Ec<i:  dynamic  compression  modulus  of  elasticity. 

Etd=  dynamic  tensile  modulus  of  elasticity. 

Etsm:  dynamic  tensile  modulus  of  elasticity  of  notch  tip 

element . 

1 

ftd:  dynamic  tensile  strength. 

Gf:  fracture  energy  density 

a:  crack  closure  stress  shape  factor.  (Fig.  3.5  and  Eq. 

3.7) 
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Figure  4.14:  Strain  Rate  Sensitivity  of  Fracture  Energy  Density 


The  crack  closure  stresses  in  the  FPZ,  which  trail  the  crack 
tip,  thus  increase  with  increasing  strain  rates.  This  result  is 
consistent  with  the  findings  for  slow  loading  rates  reported  in 
[70]  . 


As  shown  in  the  legend  for  Figure  4.13,  the  calculated 
maximum  crack  mouth  opening  displacements  (CMODmax)  the 
precast  notch  tip  never  reached  the  critical  COD,  CODc,  for 
average  strain  rates  less  than  0.0981  sec-1.  For  the  average 
strain  rates  of  0.0981  and  0.242  sec-1,  CMODmax  was  approximately 
equal  to  CODc  and  CMODmax  was  greater  than  CODc,  respectively. 
The  FPZ  was  never  fully  developed  at  average  strain  rates  equal 
or  less  than  0.174  sec-1  and  the  two  halves  of  the  fractured 
specimen  remained  rigidly  attached  to  each  other  after  the 
microcrack  had  penetrated  through  the  last  ligament  of  the 
specimen.  The  crack  closure  stress  due  to  aggregate  bridging  had 
to  be  overcome  in  order  to  separate  the  two  halves . 

Figures  4.3,  4.4,  4.7  and  4.8  show  that  the  peak  load  was 
achieved  when  the  strain  at  gage  SG03  changed  from  compressive 
to  tensile  strain.  The  numerical  results  showed  that  the 
microcrack  tip  had  then  extended  to  about  43.1mm  (1.7  in.)  from 
the  precast  notch  tip.  This  FPZ  length  of  43.1mm  (1.7  in.)  was 
independent  of  the  notch  tip  strain  rate  and  must  therefore  be  a 
function  of  specimen  geometry. 
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Since  the  FPZ  did  not  fully  develop  in  the  specimens 
subjected  to  the  slower  strain  rate  loadings,  i.e  average  strain 
rates  equal  or  less  than  0.174  sec-1,  larger  size  specimens  are 
needed  for  comprehensive  studies  of  strain  rate  effects  on  crack 
closure  stresses.  Also  slight  discrepancies  in  the  measured  and 
computed  strains  after  passage  of  the  microcrack  tip  through  the 
location  of  SG02  suggest  that  further  decrease  in  the  crack 
closure  stress  in  the  trailing  portion  of  the  constitutive 
relation  in  Figure  4.14  may  be  warranted.  Finally  the 
difficulties  with  determining  the  response  for  the  ligament  at 
the  precast  notch  tip,  due  to  shrinkage  effects,  demonstrate 
clearly  the  need  to  minimize  factors  that  cause  those  effects 
due  to  laboratory  testing  and  the  likely  importance  of  those 
effects  for  practical  applications  of  fracture  mechanics  to 
concrete  structures. 

4 . 5  Conclusions 

This  dynamic  fracture  study  of  concrete  three-point  bend 
specimens  showed  that  for  the  range  of  strain  rates  studied 
(0.0009  through  0.24  sec-1): 

(1)  The  log  of  the  tensile  strength  increased  almost  linearly 
with  the  log  of  the  strain  rate.  The  strength  at  the 
fastest  rate  of  0.24  sec-1  was  110  percent  greater  than  the 
static  tensile  strength. 


73 


(2)  The  logs  of  the  tensile  and  compressive  moduli  of 
elasticity  increased  almost  linearly  with  the  log  of  the 
strain  rate.  Tensile  and  compressive  moduli  of  elasticity 
at  the  fastest  rate  of  0.24  sec-1  were  60  and  40  percent 
larger, respectively,  than  the  values  for  static  loading. 

(3)  The  average  crack  velocity  increased  with  increasing  strain 
rate  and  reached  132  m/sec  (5200  in/sec)  for  a  strain  rate 
of  0.24  sec-1 . 

(4)  A  fully  developed  fracture  process  zone  (FPZ)  was  not 
achieved  in  any  of  the  specimens  fractured  at  strain  rates 
equal  to  or  less  than  0.174  sec-1.  In  all  specimens,  the  FPZ 
extended  about  43  mm  (1.7  in)  prior  to  the  achievement  of 
the  maximum  load  and  that  length  was  not  strain  rate 
sensitive.  Specimens  with  greater  overall  depths  than  those 
tested  here  are  necessary  for  studies  of  FPZ  length 
variations . 

(5)  The  crack  closing  stress  for  a  given  crack  width  was  strain 
rate  dependent.  However,  since  the  tensile  strength 
increased  with  strain  rate  and  the  critical  opening 
displacement,  CODc,  for  zero  closing  stress  decreased  with 
strain  rate,  the  fracture  energy  density  was  strain  rate 
independent . 
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CHAPTER  5 


IMPACT  LOADING  THREE-POINT  BEND  TESTS 
WITH  MOIRE  INTERFEROMETRY 


5 . 1  Introduction 

In  order  to  examine,  in  further  detail,  the  dynamic 
fracture  response  of  concrete  under  dynamic  loading,  one  series 
of  the  three-point  bend  specimens  reported  in  Chapter  4  was 
repeated  using  a  drop  weight  set-up  and  with  the  transient  moire 
interferometry  fringes  surrounding  the  rapidly  propagating  crack 
recorded.  This  chapter  describes  the  results  of  those  tests  and 
their  interpreation . 

The  test  set-up  was  described  prviously  in  Section  3.2 
along  with  details  of  the  Moire  Interferometry  instrumentation. 
The  experiments  used  three  point  bend  specimens  whose  geometry 
and  concrete  properties  duplicated  those  of  the  specimens 
described  in  Chapter  4. 

5 . 2  Numerical  Analysis 

The  input  conditions  for  the  dynamic  finite  element 
propagation  analysis  of  the  test  results  were  experimentally 
determined  impact  load  and  an  assumed  constitutive  relationship 
which  relates  the  crack  closure  stress  and  the  crack  opening 
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displacements,  COD's,  in  the  FPZ.  The  initial  form  used  for  this 
constitutive  relation  was  the  same  as  that  derived  in  the  prior 
displacement  controlled  tests .  Crack  extension  was  accomplished 
by  advancing  the  micro-crack  tip  i.e.  the  tip  of  the  FPZ,  when 

I 

the  crack  tip  element  reached  the  dynamic  tensile  strength,  f 

I 

for  the  concrete.  f^d  was  obtained  by  multiplying  the  static 
splitting  tensile  strength,  f  =  6.5  V  f '  c,  by  the  ratio  of  the 
dynamic  to  static  peak  loads  for  failure  of  the  SEN  three  point 

I 

bend  specimens.  The  resultant  f  ,  which  is  given  in  Table  5.1, 

I 

was  13.1  MPa  (1,900  psi) .  There  was  an  unique  value  of  f  since 
only  one  drop  weight  height  was  used. 

The  macro-crack  tip,  (the  tip  of  the  traction  free  crack 
which  is  preceded  by  the  FPZ),  was  advanced  when  the  crack 
opening  displacement  (COD)  reached  a  specified  critical  crack 
opening  displacement  (C0Dc) . 

A  continuous  function  was  used  to  relate  the  crack  closure 
stress,  fees/  and  the  COD  as  follows: 

fees  =  f"  ■  (1-C0)a  (5.1) 

td 

where  0)  =  C0D/C0Dc  is  the  crack  opening  displacement  rate,  a,  is 

the  crack  closure  stress  shape  factor  and  f  is  the  "maximum" 

td 

« 

dynamic  crack  closure  stress  which  is  equal  to  or  less  than  f 

M  td 
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TAB  LX  5 . 1 

Mechanical  properties 


*  Strength 

Static  compressive  strength;  (measured) 

=  55.2  MPa(7,990  psi) 

Static  splitting  tensile  strength;  (measured) 
ft3  =  3.96  MPa  (575  psi)  =  6.43  Vf  'c 

Dynamic  tensile  strength  (for  micro-crack  extension 
criterion)  ; 

f *  =13.1  MPa  (1, 900  psi) 


where  f  =  5.24  MPa  (760  psi),  a  -  1.5  and 
CODc  =  0.0762  mm  (0.003  in) 

*  rracture  energy  density 

Gf  =  •  f0  •  CODc  =  160  N/m  (0.912  lb/in) 

*  Crack  Velocity  (Average  value  after  12.7  rum  crack 
extension) 

a  =  290  m/s  (977  fps) 
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Appropriate  values  for  a,  CODc,  f  were  obtained,  as  explained 

td 

previously,  by  numerically  matching  the  measured  and  predicted 
global  responses  for  the  specimen.  An  iterative  analysis  was 
required  during  each  time  increment  of  the  dynamic  finite 
element  analysis  in  order  to  obtain  a  kinematically  admissible 
COD  and  a  dynamically  admissible  crack  closure  stress  along  the 
FPZ . 


The  fracture  energy  density,  Gp,  can  then  be  computed  by 

C°DC 

GF  =  J  fees  d  (COD )  =  ^  ftd  CODc  (5.2) 

M 

The  initial  guess  for  the  maximum  crack  closure  stress,  ftd' 
was  obtained  using  the  COD  variation  near  the  crack  tip.  In 
contrast  to  the  findings  in  Chapter  4,  numerical  experimentation 
showed  that  for  these  experiments  the  maximum  crack  closing 

n  * 

stress,  f  ,  could  not  be  equated  to  the  tensile  strength,  f^. 
The  physical  significance  of  this  finding  will  be  discussed 
later.  After  the  numerical  iteration  process,  the  fracture 
properties,  i.e.  maximum  crack  closure  stress,  a  and  CODc,  as 
well  as  the  tensile  and  compressive  moduli  of  elasticity,  Etd  and 
Ecdf  respectively,  were  obtained  via  an  inverse  process  in  which 
these  properties  were  optimized  to  fit  the  measured  strain 
histories  for  the  three  gage  locations  shown  in  Figure  4.1. 
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5 . 3 


Result s 


5.3.1  Experimental  Results 

Four  specimens  were  fractured  with  an  impact  load  dropped 
from  a  height  of  400  mm  (15.7  in).  As  shown  in  Figure  5.1,  the 
time  variations  of  the  impact  load  were  nearly  identical  for  the 
four  tests  except  for  an  unexplained  spike  load  in  Test  No.  3. 
Figures  5.2,  5.3  and  5.4  show  the  strain  gage  times  histories  of 
the  four  tests  for  gages  SGOl,  SG02  and  SG03,  respectively. 
Despite  the  strong  similarity  of  the  load  traces,  the  strains 
for  gages  SGOl  and  SG02  in  the  tension  region  varied 
considerably  in  the  four  tests.  Thus  "average"  result  for  those 
two  gages,  which  was  consistent  with  the  strain  histories  for 
all  four  tests,  were  used  in  the  numerical  analysis.  That 
average  history  is  shown  as  a  unbroken  line  in  Figs.  5.1  through 

5.4.  Also  shown  in  Figures  5.2,  5.3  and  5.4  are  the  time  at 
which  the  FPZ  tip  passed  a  given  strain  gage  location.  The 
strain  peaks  for  strain  gages  SG02  and  SG03  in  Figures  5.3  and 

5.4,  were  taken  as  the  time  at  which  the  FPZ  tip  passed  those 
gages.  However,  such  was  not  the  case  for  the  first  strain  gage, 
i.e.  SGOl  in  Figure  5.2,  since  previous  numerical 
experimentation  showed  that  locked-in  shrinkage  stresses, 
generated  during  the  handling  and  testing  of  the  prenotched 
specimen,  altered  the  behavior  at  that  location.  Previous 
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Figure  5.1:  Measured  and  Average  Load 


Figure  5.2:  Strain  Histories  at  SG01 


Strain  Histories  at  SG02 


numerical  experimentation  also  indicated  that  dynamic  microcrack 
tip  extension  initiated  at  0.26  msec  or  at  the  first  consistent 
change  in  strain  rate  in  the  SGOl  record.  The  continued  slight 
rise  and  the  subsequent  plateau  in  the  SGOl  strain  history 
represented  continuous  increase  in  strain  in  the  crack  tip 
element  despite  the  onset  of  micro-crack  propagation  at  about 
0.26  msec.  In  contrast,  the  passage  of  the  crack  tip  at  SG02  and 
SG03  is  characterized  by  a  precipitous  drop  in  strains  with  the 
passage  of  the  microcrack  tip.  Average  crack  velocities  were  151 
m/sec  between  strain  peaks  for  SGOl  and  SG02  and  257  m/sec 
between  strain  peaks  for  SG02  and  SG03 .  This  average  crack 
velocity  is  19%  lower  than  the  crack  velocity  320  m/sec  which 
was  obtained  from  the  moire  interferometry  data  described  next. 

Figure  5.5  shows  four  typical  sequential  frames  of  the 
moire  interferometry  fringe  patterns  associated  with  a 
propagating  crack  for  a  given  trigger  delay  time.  By  varying  for 
each  of  the  four  test  beams,  the  trigger  delay  time  of  the 
IMACON  790  from  30'  to  600  micro-seconds  in  intervals  of  100 
micro-seconds  a  sequence  of  crack  opening  profiles  with  rapid 
crack  extension  was  obtained.  Figure  5.6  shows  with  symbols 
variations  in  crack  opening  displacement  (COD)  that  could  be 
determined  from  those  moire  patterns.  The  authors'  previous 
moire  interferometry  data,  involving  stable  crack  growth  in 
similar  concrete  fracture  specimens  loaded  staticallly,  [69-71] 
showed  that  COD  variations  narrowed  toward  the  crack  tip.  In 


Figure  5.5:  Moire  Interferometry  Record  of  Specimen  No. 
Frame  Interval  10  micro-seconds 


Figure  5.6:  COD  Variations  of  a  Propagating  Crack 


contrast,  in  Fig.  5.6,  the  COD  variations  in  the  vicinity  of  the 
crack  tip  suggest  a  parabolic  crack  opening  profile.  Such  a 
crack  opening  profile,  which  is  in  agreement  with  the  results  in 
[72],  suggests  a  decreasing  influence  of  crack  bridging  and  a 
gradual  transition  to  linear  elastic  fracture  mechanics  with 
increasing  strain  rates. 

5.3.2  Numerical  Results 

The  average  impact  load  and  the  assumed  constitutive 
relation,  i.e.  Equation  (5.1),  were  used  as  input  boundary 
conditions  for  the  propagation  analysis  [66]  of  the  "average" 
result  of  the  four  tests.  The  results  and  the  experience 
obtained  from  the  prior  analysis  of  the  displacement  controlled 
tests  reported  in  Chapter  Four  greatly  reduced  the  effort  in  the 
iterative  process  necessary  to  match  the  computed  and  measured 
crack  extension  and  the  three  strain  gage  histories. 

Figure  5.7  shows  the  computed  crack  extension  histories. 
The  slow  crack  velocities  during  the  initial  phase  of  dynamic 
crack  propagation  are  believed  to  be  caused  by  the  pre-existing 
damage  region  which,  as  mentioned  previously,  is  the  resvlt  of 
the  shrinkage  effects  at  the  precast  notch  tip  during  curing  and 
subsequent  handling  of  the  specimens .  The  symbols  in  the  figure 
indicated  the  measured  crack  extension  as  recorded  by  the  strain 
gage  histories  of  SG01,  SG02  and  SG03 . 
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Figure  5.8  shows  the  measured  and  computed  average  strain 
histories  at  the  three  strain  gage  locations.  Reasonable 
agreement  between  the  computed  and  measured  strains  are  noted 
prior  to  the  arrival  of  the  microcrack  tip  of  the  FPZ  at  the 
given  strain  gage  location.  The  subsequent  decrease  in  the 
measured  strains  is  also  modelled  relatively  well,  and 
particularly  well  for  SG03,  by  the  numerical  simulation.  As 
noted  previously,  the  strains  for  SG01  are  less  than  those 
expected  for  elastic  behavior  because  the  modulus  for  the  first 
two  elements  in  front  of  the  crack  tip  was  reduced  as  described 
in  Chapter  Four. 

In  Figure  5.6,  the  solid  curves  represent  the  computed  COD 
variations  with  crack  extension.  In  this  study,  unlike  in  the 
previous  stable  crack  growth  studies,  the  computed  crack  opening 
profile  shows  a  more  blunt  crack  ^ip  that  agrees  well  with  the 
crack  tip  profile  computed  from  the  moire  interferometry 
results . 

Figure  5.9  shows  the  relation  between  the  average  measured 
load  and  the  computed  load-line  displacements.  Strain  softening 
of  the  curve  was  not  predicted  until  a  load  line  disple. '  -ment  in 
excess  of  0.25  mm  (0.01  in).  However,  since  the  crac-:  tip  was 
only  three  nodes  away  from  the  impact  point  at  that  time, 
computation  was  terminated  at  0.20  mm  '0.008  in). 
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Figure  5.9:  Measured  Load  versus  Computer  Load-Line  Displacement 


Figure  5.10  shows  the  variations  in  the  different  energies 
with  crack  propagation.  The  sum  of  the  three  energies,  i.e. 
strain,  kinetic  and  fracture,  was  at  the  worst  within  3  percent 
of  the  corresponding  input  work  during  the  dynamic  fracture 
process.  The  excellent  agreements  between  the  instantaneous 
measured  and  computed  total  energies  provides  additional 
verification  of  the  numerical  analysis. 

5.3.3  Strain  Rate  Sensitivity  of  Mechanical 
Properties 

The  final  material  properties  gleamed  from  the  entire 
numerical  analysis,  are  summarized  in  this  section.  The  computed 

I 

dynamic  tensile  strength,  ftd,  and  the  tensile  and  compressive 
moduli  of  elasticity  are  shown  in  Table  5.1.  Also  shown  are  the 
static  properties  obtained  from  compression  tests  on  152.4  x 

304.8  mm  (6  x  12  in)  cylinders.  The  increase  in  the  dynamic 

»  , 
tensile  strength,  ftd,  at  these  high  strain  rates  is  consistent 

with  the  increased  number  of  trans-aggregate  fractures  seen  on 

the  failure  surfaces  of  these  specimens  as  compared  to  the 

surface  of  the  statically  loaded  specimens.  The  latter  exhibited 

many  more  failures  at  the  aggregate-to-paste  interface.  The 

factor  causing  the  rapidly  propagating  crack  to  penetrate 

through  the  higher  strength  aggregates  rather  than  around  them 

is  still  an  unanswered  issue. 
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Computed  Energy  Partition 


Figure  5.11  shows  the  variations  in  the  constitutive 
equation,  which  relates  the  crack  closure  stress  and  the  COD  in 
the  FPZ,  as  determined  from  the  numerical  analysis.  The  crack 
closure  stress  shape  factor,  a=1.5,  is  related  to  the  slope  of 
the  constitutive  relation.  The  maximum  crack  closure  stress  of 

"  i 

f  =  760  psi  is  40%  of  the  maximum  tensile  strength  of  ftd  - 
1900  psi  necessary  to  advance  the  micro-crack  tip  of  the  FPZ. 

5 . 4  Discussion 

The  two  tier  tensile  strength  criterion  used  in  propagating  the 
FPZ  in  this  study  resembles  the  two  parameter  fracture  model  for 
concrete  proposed  by  Jenq  and  Shah  [75]  .  While  a  constant 
tensile  strength  criterion  was  used  in  this  finite  element  model 
(FEM)  to  advance  the  micro-crack  tip  of  the  FPZ,  the  same  FEM 
analysis  can  also  be  used  to  compute  the  variation  in  the  stress 
intensity  factor.  The  coarse  FEM  mesh  used  in  this  study 
required  the  use  of  a  calibrated  crack  tip  stress  field  to 
compute  the  stress  intensity  factor  variation  with  crack 
extension.  The  result  is  shown  in  Figure  5.12.  The  dynamic 
stress  intensity  factor  remained  essentially  constant  during 
crack  propagation  except  for  the  initial  and  terminal  phases. 
That  lack  of  constancy  means  that  use  of  a  critical  stress 
intensity  factor  as  a  fracture  criterion  is  not  justified  for 
these  tests.  The  tensile  strength  criterion  of  crack  extension 
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used  in  this  analysis,  on  the  other  hand,  is  enforced  apriori  in 


this  study  and  thus  f  d  remains  constant  for  all  phases  of  crack 
propagation . 

The  six  fold  increase  in  the  strain  rate  from  that  used  in 
Chapter  Four  resulted  in  the  tensile  strength  ffcd  being  equal  to 
13.1  Mpa  (1900  psi)  and  therefore  60  percent  higher  than  the 
7.93  MPa  (1150  psi)  reported  in  Chapter  Four.  That  change, 
however,  was  accompanyed  by  rapid  decreases  in  the  FPZ  strength 
and  length,  and  those  changes  while  unexpected,  are  consistent 
with  the  trend  for  the  variations  in  the  FPZ  with  increasing 
strain  rate  observed  previously. 

5 . 5  Conclusions 

The  following  conclusions  can  be  drawn  from  the  results  reported 
in  this  study. 

(1)  Variations  in  COD  data  obtained  by  moire  interferometry  can 
provide  additional  insight  on  the  dynamic  fracture 
mechanism  associated  with  a  rapidly  propagating  crack  in 
concrete . 

(2)  For  strain  rate  achieved  in  this  study,  1.48  sec-1,  the 
influence  of  FPZ  on  the  fracture  process  is  reduced 
considerably  compared  to  that  for  static  loading.  At  the 
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1.48  sec-1,  the  dynamic  fracture  process  approaches  that  of 
linear  elasticy  fracture  mechanics. 

» 

(3)  The  tensile  strength,  ffcd,  and  maximum  crack  closure  stress, 

ft 

f  .  differed  by  a  factor  of  2.5. 

td 
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CHAPTER  SIX 


CRACK-LINE,  WEDGE-LOADED  DOUBLE-CANTILEVER  BEAM  TESTS 

6 . 1  Experimental  Setup 

Crack-line  wedge-loaded  double-cantilever  beam  (CLWL-DCB) 
specimens,  with  proportions  based  on  those  described  in  ASTM  E 
561-81  [76]  for  R-curve  testing  in  metals,  were  used  to  provide 
data  for  proportions  of  more  practical  size  than  the  proportions 
of  the  three-point  bend  specimen  used  for  the  investigations 
described  in  Chapters  4  and  5.  The  geometry  for  those  CLWL-DCB 
specimen  is  shown  in  Figure  6.1.  Similar  CLWL-DCB  specimens  have 
been  used  in  static  fracture  studies  by  Jeang  [13]  and  Du  [74]. 

The  CLWL-DCB  specimens  were  selected  for  studies  of  strain 
rate  sensitivity  for  the  following  reasons: 

1.  CLWL-DCB  specimens  can  simulate  conditions  in  a  practical 
size  beam  as  shown  in  Figure  6.2.  The  distance  between  the 
center  of  the  hole  and  the  edge  in  the  CLWL-DCB  specimen  is 
about  the  same  as  the  distance  between  the  bottom  of  a  beam 
and  the  center  of  its  longitudinal  reinforcement.  The 
degree  to  which  the  wedge  is  driven  into  the  hole  in  the 
specimen  can  also  be  said  to  simulate  the  effective  degree 
of  restraint  provided  by  the  reinforcement. 
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Figure  6.2:  Correspondence  between  CLWL-DCB  Specimens  and 
Conditions  in  a  Reinforced  Concrete  Beam 


2. 


With  CLWL-DCB  specimens  a  testing  machine  can  be  used 
economically.  With  proper  lubrication  and  system  alignment, 
a  testing  machine  producing  more  than  5  times  the  maximum 
expected  test  load  will  be  adequate  if  the  included  angle 
of  the  loading  wedge  is  6  degrees. 

3.  A  segmented  split-pin  system  is  effective  for  maintaining 
the  load  line  independent  of  rotation  of  the  specimen  arms. 

The  fixtures  for  a  CLWL-DCB  test  are  more  complicated  than 
those  for  three-point  bend  tests  since  the  latter  are  completely 
constrained  by  the  three  external  loading  points  and  the  crack 
propagates  towards  the  loading  point  from  the  initial  notch  tip. 
The  fixtures  for  a  CLWL-DCB  specimen  were  designed  to  transfer 
the  wedge  load  to  the  specimen  without  any  motion  in  the  out-of¬ 
plane  direction  and  with  free  motion  in  the  loading  wedge 
direction.  Figure  6.3  shows  front  and  side  views  of  the  fixture 
used  for  the  CLWL-DCB  specimens.  The  split-pin  loading  system  is 
shown  in  Figure  6.4. 

The  fixture  and  the  split-pin  loading  system  were  designed 
for  the  following  conditions: 

1.  The  low-taper-angle  loading  wedge,  with  an  included  slope 
of  1/18  and  a  polished  finish,  provides  a  crack  opening 
velocity  of  0.05  mm/msec  (0.002  in/msec)  at  the  loading 
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Front  View 


MTS  Te«-Slotted  Toblt 
Side  View 


Figure  6.3:  Fixture  of  CLWL-DCB  Specimens 
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wedge  location  for  a  vertical  wedge  velocity  of  0.9  m/sec 
(36  in/sec.) 

2.  The  wedge  load  is  measured  by  a  dynamic  load  cell  fixed 
vertically  between  the  wedge-load  block  and  the  left-split 
pin  (Figure  6.4). 

3.  To  minimize  the  friction  force  in  the  wedge-load  direction, 
steel  rollers  are  used  between  the  base  and  the  sliding 
plates . 

4 .  The  main  rigidity  in  the  vertical  direction  is  provided  by 
the  tee-slotted  table  of  the  MTS  load  frame  to  which  the 
base  plates  are  fixed. 

5.  Rubber  pads  are  added  between  the  sliding  plates  and  the 
CLWL-DCB  specimen  to  prevent  oscillations  in  the  out-of¬ 
plane  direction. 

6.  To  prevent  rigid  body  rotations  of  the  specimen  and  to 
guide  the  direction  of  crack  growth  so  that  it  is  straight 
from  the  notch  tip,  the  four  corners  of  the  specimen  are 
fixed  to  the  sliding  plates  with  clips. 

The  two  non-contact  capacitance  gages  were  used  to  measure 

the  two  crack  opening  displacements,  2V\  and  2V2  in  Figure  6.1. 

The  load-line  displacement,  2Vp,  was  determined  by  linear 
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interpolation  from  the  two  crack-line  displacements  2Vi  and  2V2  • 
The  external  work  history  for  the  CLWL-DCB  specimen  was 
determined  by  integration  of  the  load  and  load-line  displacement 
histories.  That  history  is  being  used  to  validate  the  numerical 
simulation  of  the  test  results  since  the  external  work  history 
is  a  scalar  quantity  which  represents  the  average  responses  of 
the  entire  system.  Three  strain  gages  were  bonded  at  a  half  inch 
offset  from  the  center  line  of  the  specimen,  as  shown  in  Fig. 
6.1,  for  the  same  purpose  as  in  the  three-point  bend  test 
specimens .  The  locations  for  gages  SG01  and  SG03  were  the  same 
for  all  tests.  The  distance  of  gage  SG02  from  SG01  was  12.7  mm 
for  some  specimens  and  38.1  mm  for  others.  The  analysis 
procedure  being  used  is  the  same  as  that  applied  to  the  three- 
point  bend  test  results  for  evaluation  of  the  significance  of 
the  load,  load-line  displacement,  and  strain  histories. 

6 . 2  Numerical  Analysis 


The  numerical 

procedure  being 

used 

to  analyze 

the 

experimental  results 

is  identical  to 

that 

described  in 

the 

previous  chapter  with  the  only  difference  being  the  specimen 
geometry.  Figure  6.5  shows  the  finite  element  model  used  to 
reduce  the  experimental  data  to  the  targeted  dynamic  fracture 
and  elastic  properties  of  the  CLWL-DCB  concrete  specimens.  At 
the  time  of  this  writing,  that  analysis  was  still  in  progress 
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Unit:  mm  Thichness  50.8 


and  only  preliminary  findings  are  described  in  this  report. 
Complete  results  will  be  presented  in  an  addendum  to  this  report 
to  be  issued  shortly. 

6 . 3  Results 

6.3.1  Experimental  Results 

A  total  of  eleven  CLWL-DCB  concrete  specimens  were  tested 
at  varying  strain  rates.  The  test  conditions  for  the  eleven 
specimens  are  listed  in  Table  6.1.  The  symbols  [1]  and  [2]  in 
Column  1  of  Table  6.1  indicate  differences  in  the  location  of 
the  second  strain  gage.  For  the  symbol  [1]  the  second  strain 
gage  was  located  127  mm  (5  in)  from  the  loading  wedge's  line, 
and  for  the  symbol  [2]  the  second  strain  gage  was  located  152.4 
mm  (6  in)  from  that  line. 

Typical  load  and  strain  histories  for  the  specimen.  Test 
No.  2,  which  was  loaded  to  failure  in  about  4  msec.,  are  shown 
in  Figure  6.6.  The  loading  rate  for  Test  No.  2  was  the  fastest 
rate  in  the  whole  test  series .  The  time  at  which  the  FPZ  tip  was 
taken  as  passing  a  given  gage  position,  is  indicated  by  a  symbol 
on  the  time  history  response  for  that  gage.  The  locked-in 
shrinkage  effects  at  the  cast-in  place  notch  tip  were  much  less 
than  those  of  the  three  point  bend  specimens,  since  as  shown  in 
Fig.  6.1,  a  saw-cut  of  12.7  mm  (0.5  in)  length  was  made  at  the 
cast-in  place  notch  tip  to  relieve  those  effects. 
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Table  6.1:  CLWL-DCB  Test  Data 


m 

mm 

mm 

mm 

msm 

■ 

nVsec 

*2 

m/sec 

a 

m/sec 

H xm 

i?  SSI 

01  [2] 

12.4 

2.54 

1.63 

2.56 

3.88 

40.7 

144 

101 

2.44 

02[2] 

11.0 

2.41 

1.43 

2.18 

3.83 

16.9 

131 

95.3 

2.36 

03[1] 

10.1 

6.55 

4.50 

5.90 

14.4 

9.07 

25.1 

23.2 

1.87 

04[2] 

9.98 

11.3 

6.15 

9.80 

20.5 

10.4 

17.8 

15.9 

2.27 

05[1] 

9.34 

40.6 

19.2 

31.4 

78.8 

1.04 

4.55 

3.84 

2.21 

06[2] 

9.14 

13.8 

7.60 

13.4 

23.3 

6.63 

19.2 

14.6 

1.91 

07[2] 

8.79  |  48.5 

30.0 

50.5 

118 

1.86 

2.85 

2.62 

1.81 

08[1] 

7.51 

286 

92.0 

162 

600 

0.18 

0.49 

0.45 

1.69 

09[2] 

7.42 

108 

71.5 

110 

239 

1.00 

1.47 

1.36 

1.53 

10[2] 

7.23 

420 

172 

406 

1136 

0.163 

0.262 

0.237 

1.68 

H[2] 

5.77 

2000 

932 

- 

- 

- 

‘ 

- 

- 

Tj  is  the  time  when  a  micro-crack  tip  passes  the  depth  of  the  strain  gage  i. 
Tp  is  the  time  at  the  peak  load  Pmax- 
W3  is  the  external  work  at  T3. 
a  is  the  average  crack  velocity. 
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Figure  6.6:  Load  and  Strain  Histories  of  Test 


Figure  6.7  shows  the  COD  gage  responses.  The  load  and  2Vi 
and  2V2  are  plotted  against  2Vp  which  is  the  load-point 
displacement  linearly  interpolated  from  2V\  and  2V2.  The  2V1  and 
2V2  data  provided  additional  data,  which  together  with  the 
strain  gage  data,  were  used  to  match  the  numerical  modelling  in 
the  hybrid  experimental-numerical  analysis.  The  locations  for 
gages  2Vi  and  2V2  are  shown  in  Fig.  6.1. 

Figures  6.8  and  6.9  show  the  corresponding  data  for 
specimen  No.  07  which  was  tested  to  failure  in  about  120  msec. 
That  rate  was  about  the  average  loading  rate  for  the  whole  test 
series.  Figures  6.10  and  6.11  show  the  corresponding  results  for 
specimen  No.  10  which  was  loaded  to  failure  in  about  1.2  sec. 
That  loading  rate  was  the  slowest  loading  rate  for  which 
complete  data  were  obtained. 

Careful  observations  were  also  made  of  the  fracture 
surfaces  of  the  specimens  and  the  deviations  of  the  propagating 
crack  from  the  centerline  of  the  specimen.  Procedures  are  being 
developed  to  describe  numerically  observed  variations.  However, 
qualitatively  it  was  clear  that,  as  described  previously  in 
Chapter  4,  the  roughness  of  the  failure  surface  decreased  and 
the  degree  of  transgranular  aggregate  fracture  increased  as  the 
loading  rate  increased.  Further  it  was  apparent  that  expected 
variations  in  the  strain  and  displacement  time  histories  were 
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Figure  6.8:  Load  and  Strain  Histories  of  Test 
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Figure  6.10:  Load  and  Strain  Histories  of  Test  10 
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often  related  to  the  initiation  of  gross  deviations  in  the 
direction  of  the  propagating  crack  from  the  centerline  of  the 
specimen . 

6.3.2  Numerical  Results 

Critical  test  results  as  summarized  in  numerical  form  in 
Table  6.1  and  experimental  responses,  as  typified  by  responses 
of  Figures  6.6  through  6.11,  are  being  used  as  the  experimental 
data  for  executing  the  inverse  analysis  to  delineate  dynamic 
effects.  The  2Vi  and  2V2  data  indicate  the  approximate  crack 
opening  profile  for  the  rapidly  propagating  crack  and  the  need 
to  satisfy  that  profile  consistent  with  the  strain  and  load  data 
have  significantly  influenced  the  appropriate  choice  for  the 
constitutive  relation  of  the  FPZ.  Preliminary  analysis  favor 
constitutive  relations  of  the  form  of  Figure  5.11  rather  than 
F igure  4.13. 

6 . 4  Discussion 

As  described  above,  the  results  of  a  preliminary  analysis 
of  the  CLWL-DCB  tests  indicate  that  the  constitutive  relation 
for  a  rapidly  propagating  FPZ  is  not  only  substantially 
different  from  that  of  a  stably  growing  FPZ  but  that  the 
limiting  crack  opening  displacement  (CODc  in  Table  5.1)  also 
diminishes  in  magnitude  with  increasing  strain  rate.  It  appears 
that  the  tensile  strength  value,  ftd  Table  5.1,  is  more  strain 
rate  dependent  than  the  maximum  value  ftd*  fees  the  crack 

closing  stress.  Should  the  form  of  the  relation  of  Figure  5.11 
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be  verified  from  the  CLWL-DCB  experiments,  then  the  results  of 
Chapter  4  must  be  re-evaluated  for  consistency  with  the  results 
of  this  investigation. 

6 . 5  Conclusions 

1.  From  the  analyses  completed  to  date,  it  is  clear  that  to 
obtain  an  adequate  understanding  of  the  dynamic  fracture 
behavior  of  concrete,  specimens  with  dimensions  large 
enough  to  allow  the  full  development  of  the  FPZ,  together 
with  precise  measurements  of  the  accompanying  displacements 
and  loads  are  needed.  The  conditions  of  full  development  of 
the  process  zone,  coupled  with  loads  that  are  within  the 
capabilities  of  normal  testing  facilities,  are  difficult  to 
achieve  with  three  point  bend  specimens .  Even  special 
facilities,  such  as  the  large  drop  weight  apparatus 
reported  in  Ref.  [3],  do  not  allow  full  development  of  the 
FPZ.  Data  obtained  using  CLWL-DCB  specimens  are  highly 
desirable  since  all  necessary  conditions  can  be  achieved 
with  such  specimens,  as  demonstrated  by  the  results 
reported  here. 

2.  For  future  CLWL-DCB  tests,  it  is  desirable  that  additional 
recording  channels  be  utilized.  The  crack  opening  profile 
could  then  be  more  accurately  measured,  and  with  more 
strain  gages  discributed  along  the  centerline  of  the 
specimen,  variations  in  the  form  of  the  FPZ  and  the 
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velocity  of  the  crack  and  the  reasons  for  those  variations, 
could  be  better  identified. 

3.  For  future  CLWL-DCV  tests,  it  is  desirable  that  crack 
opening  profiles  and  displacements  for  dynamic  loading  be 
measured  using  Moire  interferometry  techniques.  As 
demonstrated  by  the  results  in  Part  I  and  in  Chapter  5  of 
this  Part  II,  the  form  of  the  crack  opening  profile  is  an 
important  criteria  for  verifying  the  form  and  magnitude  of 
the  crack  closing  pressures. 
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CHAPTER  7 


SUMMARY 

The  objectives,  described  in  Chapter  One  of  this  report, 
were  all  achieved  except  for  the  final  numerical  analysis  of  the 
CLWL-DCB  specimens.  The  results  from  that  analysis  will  be 
reported  in  a  PhD  dissertation  by  Jung  Heum  Yon  which  is 
expected  to  be  completed  by  June  1990.  Copies  of  this  PhD 
dissertation  will  be  submitted  to  the  AFOSR  contract  monitor  as 
a  supplement  to  this  report. 

Specific  results,  achieved  during  the  course  of  this 
investigation  were  as  follows: 

1.  A  displacement  controlled  dynamic  test  system  and 
instrumentation  for  dynamic  fracture  testing  of  concrete 
specimens  were  developed. 

2.  The  advantages  and  disadvantages  of  displacement  control 
versus  impact  loading  and  the  use  of  three  point  bend 
versus  CLWL-DCB  specimens  for  dynamic  fracture  testing  of 
concrete  were  identified. 

3.  A  hybrid  experimental-numerical  procedure  for  evaluating 
dynamic  fracture  properties  of  concrete  from  test  data  was 
developed.  The  ambiguity  in  this  inverse  analysis  was 
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reduced  through  the  use  of  that  procedure  combined  with 
either  the  results  of  moire  interferometry  or  dynamic  COD 
measurements . 

4.  The  strain  sensitivities  of  the  elastic  properties  and  the 
fracture  process  zone  characteristics  of  concrete  were 
determined . 

a)  The  tensile  strength  increased  with  increasing  strain  rates 
with  the  log  of  the  increase  in  strength  increasing  almost 
linearly  with  the  log  of  the  strain  rate  for  values  up  to 
about  0.2  sec-1.  For  higher  strain  rates,  the  increase  was 
greater  than  for  the  logarithmic  relation  valid  up  to  0.2 
sec-1.  At  a  strain  rate  of  0.24  sec-1,  the  tensile  strength 
was  110  percent  of  the  static  strength.  At  a  strain  rate  of 
1.48  sec-1,  the  tensile  strength  was  330  percent  of  the 
static  strength. 

b)  The  tensile  and  compressive  moduli  of  elasticity  increased 
with  strain  rate  but  percentage  increases  were  much  less 
than  percentage  increases  in  tensile  strength. 

c)  The  average  crack  velocity  increased  with  increasing  strain 
rate.  Values  reached  132  m/sec  and  230  m/sec  for  strain 
rates  of  0.24  and  1.48  sec-1  respectively.  Both  velocity 
values  are  considerably  less  than  the  theoretical  terminal 
velocity  for  crack  propagation  in  concrete. 
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d)  The  fracture  energy  density  was  found  to  be  independent  of 
strain  rate  and  equal  to  160  N/m  for  the  concrete  used  in 
these  tests. 

e)  The  influence  of  the  FPZ  on  the  fracture  process  decreased 
as  the  strain  rate  increased.  At  1.48  sec-1,  the  dynamic 
fracture  process  approached  that  for  linear  elastic 
fracture . 

f)  The  form  of  the  FPZ  changes  with  increasing  strain  rates. 
At  the  crack  tip  there  is  effectively  a  decoupling  of  the 
stress  for  advancement  of  that  tip,  the  tensile  strength, 
and  the  crack  closing  pressure  immediately  behind  that  tip. 
At  a  strain  rate  of  1.48  sec-1,  there  was  a  difference  of 
2.5  in  those  values.  The  tensile  strength  is  more  strain 
rate  dependent  than  the  maximum  crack  closing  pressure.  The 
COD  at  which  the  crack  closing  pressure  vanishes  also 
diminishes  with  increasing  strain  rate. 

g)  The  percentage  of  transaggregate  fractures  increases  as  the 
strain  rate  increases  so  that  the  effective  roughness  of 
the  fracture  surface  decreases  as  the  strain  rate 
increases . 
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CHAPTER  9 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

9 . 1  Justification 

Damage  analysis  of  concrete  based  on  dynamic  fracture 
mechanics  is  still  in  its  developmental  stage.  The  changes 
observed  in  the  fracture  process  zone  with  strain  rate  in  this 
report  vividly  demonstrate  the  decreasing  influence  of  the 
dominant  fracture  energy  dissipation  mechanism,  i.e.,  the 
fracture  process  zone,  with  increasing  dynamic  effects. 
PresumabJ  ,  a  better  knowledge  of  the  macro-  and  micro- 
mechan^s  ..s,  which  trigger  this  effect,  would  provide  the  basis 
for  designing  a  concrete  structural  component  to  retain  its 
static  fracture  process  zone  during  rapid  crack  propagation  and 
hence  be  engineered  against  dynamic  fracture.  For  example, 
trans-aggregate  failure  during  rapid  crack  propagation  could 
possibly  be  decreased  by  using  stronger  and  larger  aggregates 
together  with  stronger  interfacial  adhesion  with  the  cement.  It 
is  well  known  that  the  t-nughness  of  concrete  can  also  be 
increased  by  the  addition  of  fibers  [3,  32,  33]  and 

theoretically  it  should  be  possible  to  predict  the  response  for 
planar  damage  in  a  fiber  reinforced  concrete  by  combining  the 
response  predicted  for  the  fibers  crossing  the  damage  plane  with 
the  response  predicted  for  the  matrix  (concrete)  crossing  the 
same  plane. 
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9 . 2  Recommended  Research. 

It  is  recommended  that  a  hybrid  experimental-numerical 
method  be  used  to  delineate  the  macro-  and  micro-mechanics  of 
dynamic  fracturing  of  concrete.  The  macro-fracture  work  should 
involve  static  and  dynamic,  mode  I  and  mixed  mode  fracture 
testing  of  CLWL-DCB  specimens  and  three-  and  four-point  bend 
concrete  specimens  using  the  displacement  controlled  test  set-up 
already  developed.  All  specimens  should  have  extensive  strain 
gage  and  COD  gage  instrumentation  utilizing  the  multichannel 
recording  systems  purchased  under  the  AFOSR  Equipment  Grant  No. 
89-0123  to  the  Department  of  Civil  Engineering.  Fractography 
analyses  of  the  fracture  surfaces  and  profilometer  measurements 
of  cross  sections  normal  to  the  crack  should  be  used  to  document 
micro-mechanics  contributions  to  the  observed  macro-mechanic 
response  and  for  the  formulation  of  a  mathematical  model  to 
predict  those  contributions  from  the  physical  and  mechanical 
characteristics,  and  the  gradation,  of  the  aggregate. 

In  order  to  extend  the  maximum  rate  of  straining  to  values 
greater  than  those  which  can  be  obtained  with  the  MTS  servo- 
hydraulic  testing  machine,  tests  should  also  be  made  using 
impact  loading.  The  initial  work  should  be  done  using  small- 
scale  CLWL-DCB  specimens  impact  loaded  in  the  drop-weight 
facility  in  the  Fracture  Mechanics  Laboratory  of  the  Department 
of  Mechanical  Engineering.  That  work  should  be  complemented  with 
the  testing  of  large-scale  CLWL-DCB  specimens  impact  loaded  in 
the  drop  weight  facility  in  the  Department  of  Civil  Engineering 
at  the  University  of  British  Columbia.  These  drop-weight  tested 


127 


specimens  should  be  instrumented  so  that  with  both  strain  and 
displacement  gages,  and  with  moire  interferometry,  the  transient 
displacement  field  in  the  vicinity  of  the  running  crack  can  be 
determined . 

Linear  and  nonlinear  finite  element  analyses  should  be  used 
to  analyze  the  experimental  resilts  to  determine  the  transient 
states  displacement  of  stress  and  strain  and  hence  the  fracture 
parameters  in  those  statically  and  dynamically  loaded  specimens. 

Experimental  procedures  should  follow  those  used  previously 
for  static  loadings  with  careful  screening  of  constituent 
materials,  and  the  simultaneous  casting  and  curing  of  the 
multiple  specimens  in  a  given  test  group.  CLWL-DCB  specimens 
should  be  tested  having  wider  variations  in  maximum  aggregate 
size  and  in  aggregate  gradation  than  used  previously.  Maximum 
aggregate  size  should  range  from  3/4  to  1/4  inch.  Careful 
attention  in  the  static  tests  should  be  given  to  "creep"  effects 
so  that  strain  rate  effects  for  the  static  loading  can  be 
defined  [14].  While  the  experimental  results  should  guide  the 
direction  of  any  subsequent  test  in  a  given  series,  the  basic 
hypothesis  of  a  dual  mechanism  of  crack  bridging  at  the  FPZ  and 
brittle  fracture  based  on  LEFM  should  be  probed  analytically 
and,  if  necessary,  an  alternate  basic  hypothesis  formulated  and 
tested  analytically.  If  any  crack  branching  effects  are  observed 
experimentally,  the  significance  of  that  branching  should  be 
investigated  analytically. 
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9 . 3  Relevancy 

In  order  to  assess  the  durability  of  a  concrete  runway, 
which  is  an  essentially  unreinforced  structure  subjected  to 
repeated  impact  loadings,  it  is  highly  desirable  to  know  the 
dynamic  material,  as  well  as  dynamic  fracture,  properties  of 
concrete.  In  addition,  fundamental  knowledge  on  the  strain  rate 
dependency  of  concrete  tensile  strength  is  also  desirable  for 
assessing  the  through-wall  shear  strength  of  a  concrete  silo  or 
the  through-roof  or  wall  shear  strength  of  slabs  for  magazines 
and  other  base  facilities  which  may  be  subjected  to  either 
conventional  or  nuclear  weapon  explosions. 

In  addition  to  the  need  for  procedures  to  assess  the 
dynamic  fracture  properties  of  concrete,  there  is  an  equal  need 
for  procedures  to  assess  the  significance  of  pre-existing  cracks 
for  estimates  of  the  residual  life  of  runways,  missile  silos, 
slabs  and  walls  and  for  establishing  strategies  and  schedules 
for  maintenance,  upgrading  and  strengthening  of  existing 
structures.  The  CLWL-DCB  specimen  developed  by  the  investigators 
for  quasi-static  testing  is  ideal  for  examination  of  those 
issues  and  their  implications  for  dynamic  loading  situations. 
The  specimen  has  a  large  enough  planar  surface  that  crack 
propagation,  crack  reflection,  and  FPZ  variation  effects  can  be 
identified. 
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